NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 







COMPUTER PREDICTION OF LARGE 
REFLECTOR ANTENNA RADIATION PROPERTIES 


by 


ALAN BOTULA 


DEPARTMENT OF ELECTRICAL ENGINEERING 
NORTH CAROLINA STATE UNIVERSITY 
Raleigh, North Carolina 


July 1980 


nis research was supported by trie National Aeronautics and Space 
dir, mist rat ion through grant NSG 1888. 



ABSTRACT 


I 

I 

I 

I 

I 

I 

I 

1 

I 

I 

I 


BOTULA, ALAN BERNARD. Computer Prediction of Large Reflector Antenna 
Radiation Properties. (Under the direction of JAMES FRANK KAUFFMAN.) 

A program for calculating reflector antenna radiation patterns 
has been rewritten and extended to include a new class of reflector 
surfaces. The original program, developed jointly by North Carolina 
State University and NASA, was capable of computing patterns for 
reflectors with smooth analytic surfaces. The revised version allows 
the reflector to be composed of a number of panels, removing the 
restriction that the surface must be smooth. These Individual panels, 
however, are subject to the restrictions originally imposed upon the 
entire surface; that is, they must be analytic surfaces. 

The theoretical foundation for the program is as follows; Geomet- 
rical optics techniques are used to trace rays from a feed antenna to 
the reflector surface and back to a mathematical plane just in front 
of the reflector. The resulting tangential electric field distribu- 
tion, assumed to be the only source of forward radiation, is integrated 
numerically to calculate the radiation pattern for a desired, set of 
angles. When the reflector is composed of more than one panel, each 
panel is treated as a separate antenna, the ray-tracing procedure and 
integration being repeated for each panel. The results of the indivi- 
dual aperture plane integrations are stored and summed to yield the 
relative electric field strength over the angles of interest. 

An example and several test cases are included to demonstrate the 
use of the program and verify the new method of computation. 
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1 . INTRODUCTION 

This thesis reports the work done to extend the capabilities of 
the computer algorithm, REFLECTR, a Fortran program developed jointly 
by North Carolina State University and the National Aeronautics and 
Space Administration to compute reflector antenna radiation patterns. 

The original progiam was published In the IEEE Transactions on 
Antennas and Propagation in 1976 [1], At this stage, it was capable 
only of computing patterns for paraboloidal reflectors. It was called 
PREFPAT, an acronym for Paraboloidal REFlector PATterns. The geometrical 
formulation was later changed to provide the calculated co- and cross- 
polarized fields which would actually be measured on an antenna range. 
This revised version appeared later that year In a North Carolina State 
University report [2], Between 1976 and 1978, the program was extended 
to Include ellipsoidal and spherical surfaces. It also received its 
present name, REFLECTR. A report on these changes was published by 
NASA in May of 1978 by Dr. Pradeep K. Agrawal [3]. Dr. Agrawal con- 
tinued to improve the program, adding the ability to analyze planar 
reflectors in addition to the other surfaces [4], 

The objective of the work reported herein is to generalize 
REFLECTR to make it capable of computing patterns for more complex 
surfaces; in particular, refl ectors constructed from several panels or 
sections. The strategy for the algorithm is to represent each panel 
as a separate antenna and superpose their radiation patterns to find 
the total secondary pattern for the configuration. Among the conse- 
quences of this work were the automation of several operations formerly 
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done manually, the addition of equations for parabolic cylinders to 
allow analysis of deployable gored reflectors, a numerical technique 
for aperture boundary definition, and improved output formatting with 
more control over output operations. 

This latest version of REFLECTR Is written In standard Fortran-G 
using no compiler-specific advantages. With the exception of the 
plotting package, it should run on any Fortran compiler with little 
change. A complete description of how to organize an Input file for 
the program for several reflector configurations Is provided In 
Section 6. 

The Increase in the generality of the algorithm has not caused a 
substantial increase in computer time. In fact, the revised version 
requires only about six percent more time than the old program when 
making a pattern computation for a given single-panel reflector. 
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2. THEORETICAL DEVELOPMENT 
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2.1 Review of Theory 

The new and old versions of REFLECTR share a common theoretical 
core. Before either program Is discussed Individually, this basic 
theory will be reviewed. 

The approach taken In this algorithm Is to assume that the reflec- 
tor Is In the far field of the feed antenna and to use geometrical 
optics ray tracing to compute the electric field over a planar aperture 
just in front of the reflector surface. The tangential field over this 
aperture Is assumed to be the only source of forward radiation and is 
Integrated to yield the composite radiation pattern of the antenna. 

Edge diffraction effects are neglected. 

A reflector antenna is described In terms of two cartesian coordi- 
nate systems, the reference (unprimed) system and the feed (primed) 
system. These are shown In relation to a typical reflector antenna In 
Figure 1. The feed system is used to describe the feed antenna radia- 
tion pattern. The reflector geometry, the location of the aperture 
plane, and the overall radiation pattern are described In the reference 
system. In general, these coordinate systems are related by a rotation 
and a translation, each with three orthogonal components. 

The majority of the operations in the algorithm are performed in 
the reference coordinate system. The reflector surface is described 
by a function of the reference coordinates x, y, and z. Five types of 
surfaces are presently available to the program: planes, spheres, 

ellipsoids, paraboloids, and parabolic cylinders. Other second-order 
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polynomial functions may be added with little difficulty. The composite 
radiation pattern is expressed In a spherical coordinate system which 
Is based on the reference cartesian coordinates. Convention Is fol- 
lowed In that the theta angle is the polar angle measured from the 
positive z-axls and phi Is the angle In the xy plane Increasing from 
the positive x-axIs to the positive y-axls. 

The feed system allows the feed antenna radiation pattern to be 
described independently of Its location or orientation. The phase 
center of the feed antenna is assumed to be at the origin of Its coor- 
dinate system and the antenna points In the negative x' direction. Its 
radiation pattern is described in a spherical coordinate system, which 
Is based on the feed cartesian coordinates using the same convention as 
the reference spherical system. The feed coordinate angles are denoted 
o' and ip' to differentiate between these and the reference system 
angles. Field strength and phase information for a given set of ray 
angles is obtained by evaluating a known function, which represents an 
assumed feed pattern, or by interpolating from measured data. Because 
the feed antenna points in the negative x' direction, the main beam 
will be oriented at <p ' - 180° instead of - C ’. Hence, the radiation 
pattern must be shifted in the phi coordinate by 180° before It Is used 
in the program. 

Once a ray is generated in the feed system, it must be expressed 
In the reference system. The transformation is performed by a three- 
dimensional rotation matrix [A]. A derivation of the rotation matrix 
appears in [2], 




6 

The rotation matrix Is computed from the three angles used to 
define the rotations required to make the feed system parallel to the 
reference system. These angles will be called ALPHA, BETA, and GAMMA. 
Figure 2, taken from Agrawal [3], shows these angles as they relate to 
the feed coordinate axes. ALPHA Is a rotation about the z'-axis. It 
could be considered the azimuthal aiming angle for the feed antenna. 
BETA Is a rotation about the x'-axls, providing a means for effectively 
rotating the reflector. GAMMA is a rotation about the y'-axls. This 
angle Is commonly used to define the feed pointing angle for offset 
reflectors. All angles are defined to be positive If the rotation is 
counter-clockwise when looking In the negative direction along the axis 
of rotation. 

The reflector is Illuminated by a raster scan of rays over the 
angles it subtends. As computer storage is presently allocated, as 
many as 2750 rays, defined by their angles theta and phi, may be 
created for this purpose. This limitation may be removed by Increasing 
the dimensions of the P array (See Appendix D). For reasons explained 
later, the outermost rays are intended to miss the reflector. 

For purposes of illustration, a single ray will be traced through 
the series of operations resulting In the calculation of a point on the 
aperture plane. The ray begins as a set of angles 0 ' and in the feed 
system. The appropriate 0 '- and ((.'-polarized electric field strengths 
and the initial phase are taken from the feed antenna pattern. This 
constitutes five pieces of Information associated with a given ray. A 

A 

unit vector, s^, in the direction of the ray Is computed using the 
familiar relations 
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^ AAA 

s 1 ■ s x .x ♦ s y .y + Sj , ,2 

s x , ■ sin e' cos 

$y, ■ sin e' sin 

$ z , ■ cos 0 * 

The 0 ' and ♦' components of the field are similarly resolved to 
their rectangular components. Recall that the feed system Is trans- 
lated with respect to the reference system. The reference origin Is 
known as a point In the feed system and this coordinate triplet may be 
thought of as the translation vector (P7f expressed In the feed system 
(See Figure 1). The point where the ray strikes the reflector surface 
is defined by the vector V as shown In Figure 3. Thus, 

V • r$j - CTO 

where r Is the distance to the Intersection point from the feed origin. 
To find V, the vectors s^ and (FiJ must be expressed In the reference 
coordinate system. The transformation Into the reference system Is 
accomplished by forming a 3 by 2 matrix [BB] whose columns are the ray 

A 

unit vector s^and the translation vector, and premultiplying It by the 
rotation matrix [A]. [A] operates on each column of [BB] Independently, 

causing the ray unit vector and the translation vector to be rotated 
Into the reference system, 


[A][BB] « [B] « Cb. f j3 . 
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FIG. 3. VECTOR OPERATION 
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The resulting matrix [B] has the transformed ray unit vector as 
its first column and the transformed translation vector as Its second. 
All points in the path of the ray are now described in the reference 
system by the parametric equations 

x = b n r - b 12 

y = b 2] r - b 22 

1 = b 31 r “ b 32 

where r is the distance travelled from the feed origin along the ray. 
Simultaneous solution of the above parametric equations with the 
equation of the reflector surface will yield the point of intersection 
(xq, yg, Zq) and the distance from the feed origin R. The unit normal 
to the surface is evaluated at this point and used to compute a vector 
in the direction of the reflected ray, 

*> A AAA 

s r = s. - 2(n • s^n. 

Because the reflector is assumed to be in the far field of the feed 
antenna, the electric field incident on the reflector surface, denoted 
t. , is attenuated by a factor 1/R. The reflected field, assuming a 
perfect conductor, is given by 

A A 

E r = 2(n • E^n - £.. 

To find the intersection of the reflected ray with the aperture 
plane, another set of parametric equations is formulated, 
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x * x Q + d cos (a x ) 

y s y 0 ♦ d cos (a y ) 

2 = Z Q + d COS (a 2 ) 

where d is the distance travelled along the ray. The direction 
cosines are easily computed by 

A -4. 

(x • SJ 

cos a = 

* is r l 


cos a y 


COS a 2 


(y • s r ) 
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The parametric equations are solved simultaneously with the equa- 
tion of the aperture plane (x = x c ) to give the point of intersection 
(x fl , y a , z ) and the distance d. The phase of the field upon reaching 
the aperture plane is computed using 

♦ = MR _Ldl + * 

X 0 

where is the phase of the feed antenna radiation pattern at the 
angles defining that ray. There are five quantities computed at the 
point at which the ray intersects the aperture plane: the y and z 
coordinates, the y and z components of the electric field, and the phase 
of the field. A more detailed discussion of the above operations is 
provided by Kauffman [2]. 
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The field strength radiated by a reflector antenna at a particular 
point In space Is computed by evaluating two Integrals over the aper- 
ture plane. It Is necessary to Integrate only those points which result 
from reflections from the actual surface and not from Its mathematical 
extension. To do this, the reflector boundary must be defined on the 
aperture plane allowing the extraneous points to be Identified. Methods 
for finding and describing the boundary are discussed later. A series 
of edge points is interpolated on the boundary using information from 
points which lie outside of the aperture. The reason for over- 
illuminating the reflector as mentioned earlier is to insure that there 
will be points external to the aperture with which to interpolate edge 
points. The points external to the aperture are then discarded. The 
next processing step is to quantize the aperture plane points in their 
y-coordinate. A series of equispaced constant-y lines, called grid 
bars, is created across the bounded aperture. If the field strength 
and phase vary slowly across the aperture compared to the grid bar 
spacing, little error is incurred if the point is simply assigned a new 
y-coordinate corresponding to the nearest grid bar. 

The integrals to be evaluated are: 

E 0 =//ER 2 cos <{> exp{j[ky sin 9 sin <j> + kz cos 9 - 4>]}dydz 
Aperture 

=//(ER sin 9 + ER cos u sin 4>)exp{j[ky sin 9 sin <t> 

<p y z 

Aperture 


+ kz cos o — 4*3 } dydz. 
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The development of these Integrals may be found In Appendix A. They 
are evaluated in two steps. First, the z-dependent part of each 
integral is computed as a line integral along each individual grid bar. 
The results of the grid bar integrations, weighted by the grid spacing, 
are then summed to give the total integrals. This procedure is 
repeated for each pair of angles in the desired pattern. When the 
pattern is complete, the field value of largest magnitude in each 
polarization is assigned the value of zero decibels. The remaining 
field strength values are tabulated and plotted in dB relative to these 
maximum values. The above operations form the core of both algorithms. 

2.2 Review of the Original Program 

The latest version of REFLECTR prior to this work was organized as 
shown in Figure 4. The subroutines in the double-lined boxes will be 
briefly discussed to describe how the various tasks were distributed in 
the program and by what means they were accomplished. The main proce- 
dure (MAIN) was little more than a calling program to the five principal 
subroutines shown below it and will receive no further attention. 

NPUT was the first routine called by the main program. Its primary 
function was to read the input data for the reflector under computatioa 
In addition to the basic geometrical data, the old version required the 
programmer to supply parameters for an aperture plane boundary ellipse; 
the grid bar location and spacing; and the minimum, maximum, and incre- 
mental values of e' and <j>* for the illumination pattern. All of these 
had to be calculated prior to running the program. NPUT passed the 
illumination pattern specifications to the subroutine called FILL. 



STRUCTURE OF ORIGINAL VERSION 
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The function of FILL was to load the main working array of the 
program, hereafter referred to as the P array, with illumination data. 
FILL had no established form. It was simply a name assigned to any 
subroutine which would accept minimum, maximum, and Incremental values 
of theta and phi and load the P array with ray data. If some type of 
Information was needed for this purpose, such as a table of field 
strengths for the feed antenna, it was read in at this time. FILL 
would load the P array and return to NPUT, which would then compute 
sane trigonometric constants and return control to the main procedure. 

The next subroutine called was APERTUR. APERTUR was responsible 
for tracing the rays to the aperture plane, determining whether or not 
they resided within the aperture boundary ellipse, and interpolating a 
set of edge points. The routine began by computing the rotation 
matrix, then proceeded to trace each ray to the aperture plane. When 
the location of a point on the aperture plane was found, it was tested 
for residence by evaluating the aperture plane ellipse equation using 
Its y- and z-coordinates. Aperture blockage was tested similarly. If 
a point was inside the ellipse, it was stored at the beginning of the 
P array. If a point was on the edge of the ellipse, it was stored at 
the end of the array. If a point was found to be outside, it was dis- 
carded. When two successive points were on opposite sides of the 
ellipse, an edge point was interpolated between them and stored with 
the naturally occuring edge points. 

QUANTIZ prepared the aperture plane for integration by quantizing 
the locations of the aperture plane points in the y-coordinate. The 
edge points found by APERTUR were used in a spline interpolation 
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routine to compute a new set of edge points at the grid bar locations. 
The interior points were quantized separately and concatenated in the 
array with the new edge points. Finally, the array was sorted; first 
with respect to y, then with respect to z. 

APRPLT made a plot of the locations of the computed points on the 
aperture plane both before and after quantizing. Different symbols 
were used for every three dB variation in electric field strength 
across the aperture. A modified version of APRPLT is used in the 
revised program. 

A feature of the old program was that in the radiation pattern 
request, one could specify a pattern referenced to the beam maximum. 

The computer estimated the beam maximum by finding the ray correspond- 
ing to 9 = 90°; <j> = 180° in the feed system and stored the angles of 
its reflection from the reflector surface in the reference system. 

This feature was not retained in the revised program. 

2.3 Numerical Techniques for the Revised Program 

In order to accommodate multi-panelled surfaces, several new 
numerical techniques were added to the algorithm. These additions are 
now discussed in detail. 

Each panel of the reflector is described by a surface code, sur- 
face parameters, and a list of edge points. The surface code and 
parameters allow the computer to select the correct surface equations 
with which to compute the ray reflections. The edge points, hereafter 
referred to as perimeter points, are the vertices of a convex polygon 
which approximates the boundary of the panel in question. The perimeter 
points are expressed in the reflector cartesian coordinate system and 


wist be Input to the computer in order as one follows the boundary of 
the reflector panel either clockwise or counterclockwise. These need 
not lie in a single plane as they will be projected onto the aperture 
plane in order to estimate the aperture boundary. 

The large number of panels possible in a reflector antenna moti- 
vated the development of an algorithm to delineate the angular extent 
and incrementation of the feed pattern to be used for the illumination 
of a given panel. By performing the reverse of the ray-tracing proce- 
dure, the perimeter points can be expressed as pairs of angles in the 
feed spherical coordinate system. The details of this operation may 
be found in Appendix B. With this information, the angles subtended by 
the panel may be found and the illumination array created. The angles 
defined by each perimeter point are stored in the illumination pattern 
with the generated angles and projected to the aperture plane. These 
projected perimeter ponts define the aperture boundary for the particu- 
lar panel . 

The illumination pattern is composed under three constraints. The 
first constraint is that a certain number of rays are to be generated. 
This number, NPTPPl, is supplied by the programmer for each panel. It 
is the approximate number of rays which he desires to use in illuminat- 
ing the particular panel. The algorithm will generate a pattern with 
equal theta and phi increments which contains a number of rays as close 
to NPTPPl as possible. Second, the algorithm will generate an addi- 
tional set of rays one angular increment beyond the angles actually 
subtended by the reflector in order to slightly over-il luminate the 
panel. These extra rays are reflected from the mathematical extension 
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of the surface and used to Interpolate edge point data later In the 
program. The third constraint Involves quantizing considerations. If, 
In the feed system, phi Is held constant and theta Is varied, a 
constant-phi contour will be traced on the aperture plane. The program 
attempts to make the grid bars fall on these constant-phi contours. 

When this is done, the distance the points must be moved during quantiz- 
ing Is minimized. Another consideration concerning the grid bars Is 
that the outermost bars come close to the edge of the reflector. If 
this is not done, the panel will be represented In quantized form as 
being narrower than is actually the case. If both of these conditions 
are to be met, some forethought must be exercised when creating the 
illumination array. The illumination pattern is constructed such that 
the outermost phi angles exceed the angular extents of the reflector 
panel (as determined by the perimeter points) by 4/5 of an angular 
increment. The angles just within the angles subtended by the reflec- 
tor are then 1/5 of an angular increment inside the boundary angles. 

The resulting point distribution on the aperture plane will then 
exhibit qualities favorable to quantizing. The creation of the 
illumination pattern is thus a more complex procedure than it was In 
the old program, but it is almost completely automatic and less haphaz- 
ard. It should be noted here that this system is hardly infallible 
and tends toward greater deviation as the rays become less parallel 
with the x-axis of the reference system; that Is, for panels which lie 
near the edges of large reflectors. In these regions, the constant-phi 
contours bend into arcs, and a grid bar may cross several contour line* 
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The field strength Is usually low In these regions, however, and this 
deviation has not been found to cause problems. 

A major difference exists between the old and revised programs 
regarding the method by which the aperture plane points are prepared 
for Integration, both In the handling of the edge points and In the 
determination of the aperture plane boundary. No attempt Is made to 
distinguish between Interior and exterior points In the APRTUR routine; 
Instead, all of the rays are traced to the aperture plane and the array 
of points Is delivered to the quantizing routine. The quantizing 
routine computes the position and spacing of the grid bars In accor- 
dance with the criteria stated earlier and assigns each point, except 
for those which are the projections of the perimeter points, a new 
y-coordlnate corresponding to that of the nearest grid bar. 

The aperture boundary is convex and closed; thus, each grid bar 
will Intersect It at two points. The Intersection points may be con- 
sidered the endpoints of the grid bar, since any aperture plane point 
not lying on the resulting line segment does not reside In the aperture 
Once the z-coordlnate of each endpoint is known, residence of a data 
point In the aperture may be determined by comparing the z-coordlnate 
of the data point with those of the endpoints for that particular grid 
bar. 

To find the endpoints, the computer works Its way around the polyg- 
onal boundary, formulating the line equations given by adjacent 
perimeter points. When the line equation for a particular pair of 
points has been found, any grid bars lying between the generating points 
are intersected with the line and are assigned the resulting z-coordlnate 




20 


as one of their endpoints . After the entire boundary has been 
processed, each grid bar will have two endpoints and the extraneous 
points may be Identified. This operation Is pictured In Figure 5. 

One of the major changes In the program Is the procedure of Inte- 
grating the aperture field panel by panel Instead of all at once. 

T h1s change Is numerically justified by the fact that the phase of the 
fields for various panels Is measured In absolute angle from the same 
phase reference and the Integrations are all performed over the same 
plane. It Is justified ^rom a theoretical standpoint In that the 
linearity of the Integra's and Maxwell's equations allow the superposi- 
tion principle to be applied. 

The program has the ability to compute single-panel antennas using 
an elliptical rather than a polygonal boundary. This has the advantage 
that a smooth edge will be created and only four perimeter points are 
required In the Input data. The points supplied to the program are 
those constituting the maximum and minimum y- and z-coordl nates; that 
Is, the highest, lowest, leftmost, and rightmost points on the reflec- 
tor as seen by the feed antenna. These points are projected to the 
aperture plane in the same way as those used to define a polygonal edge, 
and an ellipse Is formulated using the distance between the high and 
1o>' ^-coordinate' as the minor axis and the distance between the high 
and low y-coordlnates as the major axis (See Figure 6). The grid bar 
endpoints are computed by solving the ellipse equation for z using the 
y-values of the grid bars. The solution Is multi-valued, and the two 
solutions for each grid bar are its endpoints. This option will cause 
the edge to be misrepresented in configurations where the Illuminating 
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rays take on extreme values of the angular coordinates, but It does 
serve as a valuable simplification In the use of the program In many 
Instances. 
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3. THE STRUCTURE OF REFLECTR 


3.1 Overview 

A block diagram of the revised version of REFLECTR Is shown In 
Figure 7. The primary subroutines, pictured In the double-lined boxes, 
will be examined In detail In this section. The main procedure Is last 
since a knowledge of the other primary routines is prerequisite to Its 
discussion. The remaining subroutines function mainly in a support 
capacity and will be discussed In Section 5. 

With the exception of the plotting package (APRPLT and PLTSET), 
the entire program Is written using double-precision constants and 
functions. IBM double precision employs eight bytes of storage for 
each real variable and sixteen bytes for each complex variable. This 
corresponds to fifteen significant figures in real variables and In 
both the real and Imaginary parts of the complex variables. The pro- 
gram has been tried with four-byte variables with a significant degra- 
dation in the results. In its present form, REFLECTR requires 258K of 
memory . 

As a means of introducing some of the variables in the program, 
the eight common storage blocks will be listed and explained: 

BLOCKG/YCBL, ZCBL, HFMABL , HFMIBL 

(Aperture plane blockage information.) 

YCBL, ZCBL - 

y and z center coordinates of the aperture plane blockage 
ellipse. 

















HFMABL, HFMIBL - 

Half-major and half-minor axes of the aperture plane blockage 
ellipse (along y and z, respectively). 

C0NTRL/N0PT(3) , NLIST, IOPT, ItlST(lOO) 

(Program control variables.) 

NOPT ( 3 ) - 

Three numbers specifying options regarding printer, plotter, 
and aperture plane data output, respectively (See Section 4). 

NLIST - 

The number of panels for which the algorithm will print 
complete Illumination and quantizing data. 

IOPT - 

A variable which Is zero when the program Is to run normally 
and one when the single-panel option Is in effect. 

ILIST(IOO) - 

The specific panels for which the algorithm is to provide 
complete Illumination and quantizing data (See Section 4). 

DIMENS/YDIM, ZDIM, YCT, ZCT 

(Dimensional data for the bounded aperture.) 

YDIM, ZDIM - 

The maximum y- and z-dlmensions for the aperture to be 
plotted by APRPLT or APRMAP. 

YCT, ZCT - 

The y and z center coordinates of the aperture to be plotted. 

EXTENT/YMIN, YMAX, ZMIN, ZMAX 

(Dimensional data for the bounded aperture. ) 

YMIN, YMAX - 

The lower and upper limits of the Lounded aperture in the 
y-coordinate. 
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ZMAX , ZMIN - 

The lower and upper limits of the bounded aperture In the 
z-coordinate. 

(This Is actually the same Information as Is contained In the 
DIMENS block but In a form which Is more convenient for 
operations other than plotting.) 


FEED/EP(91 ) , ET(91 ) , NP, NT. 

( Feed antenna parameters . ) 

EP(91 ) , ET(91 ) - 

Arrays containing the electric field strengths of the feed 
antenna In one degree increments off-axis in the e * 90° 
and 41 * 180° planes, respectively. 

NP, NT - 

The number of increments of phi and theta used In the illu- 
mination pattern, respectively. 

(The contents of this block will change depending on how the 
feed antenna is described. The above block was used In the 
first test case for which the feed pattern was interpolated 
from E- and H-plane tables.) 


MATH/PI, PI2, PID2, DTOR, RTOD 
(Mathematical constants.) 

PI, PI2, PID2 - 

Pi, two times pi, and pi divided by two. 

DTOR, RTOD - 

Constants for degree-to-radian conversion and radian-to 
degree conversion. 


PARAMS/AORORF, BELLP, PSI, ?LNPNT(3) , PLN0RM(3), FEED(3) , ALPHA, BETA, 
GAMMA, XL AM, SURFCE, NPNL, NPOINT 


(Antenna system parameters.) 


AORORF - 


The focal length of a paraboloidal reflector » the focal length 
of a parabolic cylindrical reflector, the radius of a 
spherical reflector, or the semi -major axis of an ellipsoidal 
reflector (along x), depending on which surface Is Intended 
to represent the reflector. 

BELLP - 

The semi-minor axes (along y and 2 ) of an ellipsoidal reflec- 
tor surface. Note that this does not define a completely 
arbitrary ellipsoid since the axes along y and z must be 
equal . 


PSI - 


The angle of rotation of a parabolic cylinder whose focal 
line Is In the yz plane as measured from the positive y-axls 
when looking In the negatlve-x direction In the reflector 
coordinate system (See Appendix C). 

PLNPNT(3) - 

The coordinates of a point on a planar reflector surface 
(x, y, z). 

PLN0RM(3) - 

The components of a unit normal vector to a planar reflector 
surface (x, y, z). 

FEED(3) - 

The reference coordinate system origin as expressed In the 
feed coordinate system (x' t y' t z'). 

ALPHA, BETA, GAMMA - 

The angles of rotation between the two coordinate systems as 
described in Section 2.1. 

XLAM - 

The wavelength of the radiation emitted by the feed antenna. 
SURFCE - 

An Integer variable containing the value of the surface code. 

1 - surface is a plane. 

2 - surface Is an ellipsoid. 

3 - surface Is a sphere. 
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4 - surface Is a paraboloid. 

5 - surface Is a parabolic cylinder. 

NPNL - 

The number of panels comprising the reflector. If the value 
supplied to the computer Is one or greater, the program 
expects a list of perimeter points and other surface para- 
meters for each panel and the aperture boundary is approxi- 
mated by a polygon. Entering the value of zero causes the 
program to use an ellipse to represent the boundary— the 
single-panel option. 

NPOINT - 

The number of points or rays which are being stored for pro- 
cessing In the P array at any given time. 


PATTRN/ ET0T(2,400) , AMINOR (3,5), AMAJ0R(5), MIN0R(5), MAJ0R(5), 
NANGLE(5) 

(Information concerning the overall radiation pattern.) 

ET0T(2,400) - 

A complex array which contains the summation of the aperture 
plane integration results from the individual reflector 
panels. The first subscript is one or two, denoting the e- 
or <t>-directed electric field, respectively. The pattern cuts 
are stacked end to end in the array. Presently, 400 angles 
may be stored. 

AMIN0R(3,5) - 

A pattern request array. The first subscript Is one, two, 
or three, denoting the lower, upper, or incremental values 
of the minor angle of the pattern. The second subscript 
denotes the particular request. 

AMAJ0R(5) - 

A pattern request array which contains the angle defining 
the principle plane of the pattern cut. The subscript 
denotes the particular request. 

MIN0R(5) - 

A pattern request array containing the angular coordinate 
(theta or phi) to which AMINOR applies. The subscript 
denotes the particular request. 


30 


MAJ0R(5) - 

A pattern request array containing the angular coordinate 
(theta or phi) of the principle plane of the pattern cut. 

The subscript denotes the particular request. 

NANGLE(5) - 

Contains the total number of angles In each requested pattern 
cut as determined by AMINOR. The subscript denotes the 
particular request. 

The remaining variables will be Introduced as the subroutines are 
examined. The reader Is referred to Appendix D for a listing of the 
code while reading this section. 

3.2 NPUT 

NPUT Is primarily an Input/output routine. As It begins, IOPT, 
the single-panel option variable. Is set to zero, Indicating that the 
program Is to run In the multi-panelled reflector option. The first 
four cards are read Into the TITLE array as character data. NPUT then 
reads values for FEED, ALPHA, BETA, GAMMA, XLAM, SURFCE, NPNL, AORORF, 
BELLP, PSI, PLNPNT, PLNORM , XC, YCBL, ZCBL, HFMABL, AND HFMIBL. These 
numbers are all read regardless of whether or not they are applicable 
to the reflector under study. If a value Is not needed, It Is simply 
not used. Following this, any data required for the feed antenna 
pattern is read by appropriate statements. This section of NPUT Is 
often changed to accommodate differences in the Information required by 
a particular FILL routine. Any data read here is passed to FILL through 
the FEED common block. If NPNL was entered as zero, the four extreme 
points on the reflector are read into the P array. This is followed by 
the output option data. 


Statements 77 to 98 are the pattern request Input loop. AMINOR, 
AMAJOR, MINOR, and MAJOR are read for consecutive requests during which 
the total number of pattern cuts and angles are monitored. If the 
number of cuts Is greater than five. If the total number of angles 
Involved Is greater than 400, or If any one pattern contains more than 
75 angles, an error message Is printed and the program terminates. 

These maximum numbers are storage constraints which may be changed, If 
desired, and If sufficient memory Is available In the computer. The 
number of pattern cuts Is limited by the dimensioning of AMINOR, AMAJOR 
MAJOR, MINOR, and NAN6LE. The total number of angles In all patterns 
Is limited by the dimensioning of ETOT, FIELDY, and FIELDZ. The number 
of angles In any given pattern is restricted by the dimensions of the 
arrays YFLD, ZFLD, and PWER in the main procedure. 

The DO-loop to statement 108 clears as many locations In the ETOT 
array as are required to store all of the angles in the patterns. Fol- 
lowing this step, the Information received by the computer Is printed 
and the mathematical constants in the MATH common block are calcu 1 ated. 

3.3 APRTUR and APRIN 

In the first call to APRTUR, the rotation matrix is computed from 
the rotation angles ALPHA, BETA, and GAMMA. Because it is a rotation 
matrix, [A] belongs to the class of Orthogonal Matrices. An Orthogonal 
matrix is non-singular and its inverse is equal to its transpose. AINV 
is thus obtained simply. The NPTPPL variable is set to 2250 to signify 
that in the case of the single-panel option, 2250 rays are to be gener- 
ated to illuminate the panel. NPERIM is set to four in order to accom- 
modate the four extreme points which are entered when the single-panel 
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option 1$ used. If the program Is to run normally; that Is, using a 
polygonal aperture boundary, the APRIN routine will be called and 
NPERIM and NPTPPL will be updated to the appropriate values. In sub- 
sequent calls to APRTUR, the above steps are skipped. 

If the single-panel option Is not In effect, the APRIN subroutine 
Is called to read data for the reflector panel. APRIN reads NPERIM, 
SURFCE, and NPTPPL on the first card of the panel data and screens 
these values for easily detected errors. If any are found, an error 
message Is printed and the program stops execution. The perimeter 
points are then read Into the P array with the x-, y-, and z-coordlnates 
occupying the first, second, and third locations, respectively. These 
are followed by any other geometrical information needed for that 
surface. For example, were the panel an elliptical surface, values 
would be expected for AORORF and BELLP. If the output control variables 
dictate that data Is to be printed for this panel, APRIN will do so. 

If not. It returns to APRTUR. After a call to APRIN, the values of 
SURFCE and whatever other parameters are required for the panel have 
been updated. For example. If the overall reflector was not an 
ellipsoid, an ellipsoidal panel would cause AORORF, BELLP, and SURFCE 
to be changed. 

The statements In APRTUR between 50 and 65 are concerned with 
finding the angles subtended by the reflector or the reflector panel. 
The perimeter points are expressed as angles in the feed system and 
monitored to record the maximum and minimum 0' and <{>' angles repre- 
sented by the set of points. This procedure is applied to both the 
perimeter points and the four extreme points in the single-panel 
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option. Appendix B details the geometrical formulation for this 
operation. 

The statements between 65 and 95 create the Illumination array 
from the Information found above. DELP and DELT are the angular 
dimensions of the panel In phi and theta* respectively. It Is required 
that 


NP * NT « NPTPPL. 

An estimate of an appropriate NP Is found by 

NP*[NP* ( DELT /DELP ) ] » NPTPPL 

NP - [(DELP/DELT)*NPTPPL]**0.5 + 1. 

One Is added to round NP to the next larger Integer. The statement 
following this Insures that the Illumination pattern will have a set of 
angles generated at the phi angular center of the reflector by forcing 
NP to be odd. 

If the reflector subtends DELP radians. It Is desirable to have 

ANGINC*(NP-1 ) • DELP + (0/8*ANGINC)*2. 

(See Section 2.3 concerning the Illumination pattern.) Simplifying the 
above, the angular Increment Is found to be 

ANGINC * DELP/ (NP - 2.6). 

NT is computed such that a set of rays exist at the e' angular center 
of the reflector and the rays corresponding to the upper and lower 
limits of o' miss the reflector. The resulting maxima, minima, and the 
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total number of angles generated are computed and. If desired, printed. 
The DO-loop to statement 95 loads the complete set of Illumination 
angles Into the P array just after the angle pairs corresponding to the 
perimeter points. The FILL routine Is then called to supply the angle 
pairs In the P array with field strength and phase Information. 

3.4 FILL 

Many FILL subroutines have been written for REFLECTR. One of the 
most useful Is described In this section. It may be used for any 
antenna for which the E- and H-plane patterns are available as long as 
these patterns are symmetrical with respect to the negative x'-axls and 
may be adequately represented when sampled at one-degree Intervals. 

The subroutine must be provided with two tables of relative field 
strengths In one-degree Increments— one In the e' ■ 90® plane and one 
In the $' » 180® plane. 

Consider the feed coordinate system and a typical ray as shown In 
Figure 8. The ray and the x'-axls define a plane. Let 6 be the angle 
formed between the ray and the negative x'-axls In this plane. 6 Is 
given by 


6 » arccos(s1n o' cos <j>'). 

Because It Is the angle between the ray and the axis of radiation 
symmetry, 6 determines the values for the electric field to be taken 
from both the e' * 90® and <p ' s 180® patterns. Interpolation from the 
field strength tables Is linear. Let c be the angle between the x'y' 
plane and the plane containing the ray and the x'-axls. Measured from 
the fourth quandrant of the x'y' plane, e Is given by 
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This angle determines th* proportion of each pattern used In computing 
the final electric field strength for the ray. The Interpolation 
scheme Is based on an ellipse as shown In Figure 9. The half-major 
and half-minor axes are the 6' » 90* plane and a' • 180° plane field 
strengths, respectively. The ellipse equation Is written In polar form 
to express the radius In terms of the angle c. 

( r - S — )* * ( r -* — > 2 ■ 1 
'«90 V“180 


u ■ r cos c 
v ■ r sin e 


r » e ■ -* v t. : », ™ __ 

t0t £,-.90 ‘ tE 4-.180 C ° 5t 

The code for this subroutine may be found In Appendix C. PROJX Is 
the direction cosine of the ray with the negative x'-axls; that Is, the 
cosine of the angle 6. PROJY, the cosine of the angle c, Is set to 
zero, t Is undefined when 6 = 0; thus. It Is not computed In this ca..: 
ANGIX Is the measure of 6 In degrees. L0 and IHI are the Indices of 
the array locations for the Integral angles just below and just above 
6, respectively. These are used to linearly Interpolate the electric 
field In each principal plane. The elliptical Interpolation Is then 
carried out and the result Is loaded Into the proper location - ,ie 
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P array. The subroutine Is presently written so that the antenna Is 
♦'-polarized and the Initial phase is zero. 

3.5 QUANTZ 

It is assumed that the rays representing the extreme values of ♦' 
miss the reflector; hence, there are NP-2 constant-^' contours on the 
reflector. NBARS, the number of grid bars, is set to this value. An 
appropriate number of locations in the array of grid bar endpoints, Z, 
are set to 10^. This number is used as a flag later in the program 
to determine whether the particular location has received an endpoint. 
After initializing the Z array, a loop finds the maximum and minimum 
y- and z-coordi nates represented in all of the perimeter point projec- 
tions and writes these points into the PR5 array. This array must be 
5 by n where n is one greater than the number of perimeter points 
expected for any particular panel. The array is presently 5 by 41. 

The information found in this loop is used to compute plotting para- 
meters (YCT, ZCT, YDIM, and ZDIM) and the grid bar locations. Note 
that the grid bars are placed in such a way that the outermost grid 
bars are 1/5 of a grid spacing inside the aperture. Lines 31 through 
41 constitute the quantizing loop. A point is first examined to deter- 
mine whether the y-coordinate lies within the established y-range of 
the aperture. If it does, it is assigned the y-coordinate of the 
nearest grid bar. If not, it is discarded. NPOINT is adjusted to 
account for the points which are deleted in this operation. The 
remaining points are sorted; first with respect to y, then with respect 
to z. The sorting process may be thought of as arranging the grid bars 
end to end in the array. 
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The next processing step Is to compute the grid bar endpoints. 

The method using perimeter points Is described first. As stated before, 
the basic strategy is to formulate line equations with adjacent perim- 
eter points and intersect these with any grid bars which might lie 
between them. The variables Involved are defined below: 

KDEX - 

The index of the perimeter point currently under considera- 
tion. 

Yl, Z1 - 

The y- and z-coordinates of the first point used to formulate 
the line equation. 

Y2, Z2 - 

The y- and z-coordinates of the second point used to formu- 
late the line equation. 

SLOPE - 

The slope of the line. 

B - 

The z-intercept of the line. 


YLO, YHI - 

The lesser and greater values of the y-coordinate of the two 
perimeter points under consideration. 

Z(2,71) - 


f The array of grid bar endpoints. Z(1,J) and Z(2,J) are the 

A lower and upper endpoints of the Jth grid bar, respectively. 


[ 

[ 

I 
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INDEX - 

The index of the location in the Z array in which the 
computed endpoint is stored. 
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YQ - 

The y-coordinate of the grid bar which Is to be Intersected 
with the line. 

ZEE - 

The z-coordlnate of the computed endpoint. 

I LOAD - 

A loading Index controlling which of two possible locations 
in the Z array will receive the computed endpoint once the 
second subscript is fixed. 

A flowchart of the process used to load the Z array Is shown In 
Figure 10. To begin, the first point in the PR5 array is copied Into 
the location following the last perimeter point. KDEX is set to 2 to 
start the loop at the second point in the array. Y2 and Z2 are 
assigned the y- and z-coordinates of the first point. Y1 and Z1 are 
taken from the second point in the array and the difference between Y1 
and Y2 is checked to insure that the line will not have infinite slope. 
Since the boundary is assumed to be convex, a vertical line must have 
one of the two outermost y-coordinates of the reflector. The grid bars 
were placed inside of these extremes; thus, a vertical line will inter- 
sect no grid bars. If Yl and Y2 differ by too small an amount, the 
algorithm rejects the line and goes to the next point pair. Otherwise, 
the slope and z intercept are computed and YLO and YHI are determined. 
The first grid bar beyond YLO is the first candidate for intersection 
with the line. The y-coordinate and its index in the Z array are com- 
puted. If the grid bar Is beyond YHI, the line is abandoned and the 
next line is formulated. If it is within the range of the line, the 
intersection is found and stored in the Z array. The index and 
y-coordinate are incremented by 1 and GRID, respectively, until YHI has 



FIG. 10. FLOWCHART FOR LOADING Z ARRAY 
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been exceeded. The next line Is then formulated and the process con- 
tinues until the boundary has been closed. 

After the Z array Is loaded, a short loop following the loading 
procedure Insures that Z(1 ,J) Is less than Z(2,J) and that all locations 
of the array have been loaded with intersection coordinates. 

If the single-panel option is in effect, the endpoints are found 
by solving the boundary ellipse equation at the grid bar locations. 
HFMAEX and HFMIEX are the half-major and half-minor axes of this ellipse, 
respectively. 

The next section of QUANTZ imposes the boundary on the aperture. 

The approach taken is to step through the data points in each grid bar 
in ascending order, discarding those which are outside the aperture and 
retaining those which are inside or on the edge. When a grid bar end- 
point is crossed, a new data point is interpolated with that z- 
coordinate. Any point found to reside in the aperture plane blockage 
ellipse is assigned zero field strength and phase. Upon completion of 
a grid bar, the number of retained and interpolated points is checked 
and the grid bar is deleted if this number is less than three. This is 
a requirement of the integration routine. A flowchart for this section 
of QUANTZ appears in Figure 11. The variables are defined below in 
order of their appearance. 

N - 

One greater than the total number of points retained in the 
aperture (including blocked points) at any given time. 

L - 

The total of points retained on the present grid bar 
(including blocked points) at any given time during the 
processing of that grid bar. 
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PBLK - 

An array representing a blocked point. The field and phase 
values are zero and the first two locations contain the y- 
and z-coordlnates of the point under consideration. 


YQ - 


The y-coordlnate of the present grid bar. 

IDEX - 

The second subscript of the location of the Z array which 
contains the endpoints of the grid bar at YQ. 

TEST - 

A variable which is negative, zero, or positive if the data 
point Is outside, on, or between the grid bar endpoints, 
respectively. 

TESTBL - 

A variable which is negative, zero, or positive if the data 
point is outside, on, or inside the blockage ellipse, 
respectively. 

TESTO - 

The value of TEST for the previous data point. 

POLD - 

The previous data point. 

NCHG - 

A variable controlling where a new point is loaded into the 
P array. NCHG = 0 means that the new point is to be loaded 
normally, while NCHG = 1 means that the point is to be 
loaded one location before the location indicated by N. 


The only complication In the process occurs when the algorithm 
encounters an interior point after an exterior point. A point is 
Interpolated on the boundary after the interior point has already been 
stored in the array. The interior point should follow the interpolated 
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point In the array because of Its higher z-coordlnate; thus, the points 
would be out of order If the Interpolated point were loaded normally. 
When this situation arises, the Interior point Is moved up one space In 
the array and the Interpolated point Is loaded into the gap by subtract* 
ing NCH6 from the loading Index N. If the Interpolation Is a result of 
moving from an interior point to an exterior point, NCH6 Is zero and 
the Interpolated point Is loaded In normal order. 

The end of QUANTZ sets NPOINT to the proper value and returns the 
perimeter points to the P array, storing them after the aperture points. 
Quantizing data Is then printed and the P array Is written to a disk 
file (excluding the perimeter points at the end) If the output control 
code Indicates. 

3.6 tntgrT 

INTGRT is essentially the section of the old RADPAT subroutine 
which computed the radiation pattern integrals. The sections which 
performed other tasks have been deleted or moved to other parts of the 
program. INTGRT accepts the quantized P array and the pattern request 
from the main procedure and returns the integrated results for each 
angle of the pattern in the FIELDY and FIELDZ arrays in complex form. 
The integrals evaluated by INTGRT are: 

E q = // ER Z cos 9 exp { j [ky sin e sin 4> + kz cos e - 4>]}dydz 
Aperture 

E =//(ER sin 0 + ER z cos 0 sin <f>)exp{j[ky sin e sin $ 

▼ y 

Aperture 

+ kz cos 0 - 4 >]}dydz 
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where ER y , ER 2 , end ♦ are the y- and z-components and phase of the 
electric field at the point (y,z) In the aperture. The derivation of 
these Integrals appears In Appendix D. Each of these Integrals Is 
divided Into a set of Integrations along each grid bar and a transverse 
Integration of the results from the grid bars. The grid bar Integrals 
are: 


Z! * / ER Z cos $ exp{j[kz cos 0 - *] }dz 
y * const. 

YI * / (ERy sin e + ER 2 cos 0 sin <|>)exp{j[kz cos 0 -e]}dz. 
y * const. 

The total Integrals are the transverse Integrals. 

E 0 * FLDZ * / ZI exp{jky sin e sin <j>}dy 
All bars 

* FLDY = / YI expfjky sin e sin <j>}dy 
All bars 

Integration is by the trapezoidal rule with the accuracy being enhanced 
by linearly interpolating a series of equidistant intermediate points 
between the existing data points. The number of subdivisions, called 
NPARTS, is presently seven. 

A partial flowchart for INTGRT is shown In Figure 12. The variable 
SEN signifies the beginning and end of the P array. The first time It 
Is encountered Indicates that the first grid bar has just been Inte- 
grated. The second time indicates that the last grid bar has been 
Integrated. The subroutine begins by determining the first e and <p 
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angles In the pattern and computes their sines and cosines as 
required In the Integrals. Although the Integrals are evaluated 
simultaneously, only the E 0 Integral will be discussed for the sake 
of clarity. The first two grid bars are Integrated and these results 
form the first trapezoid of the transverse Integral. This trapezoid 
Is evaluated and stored In FLDZ. The next grid bar closes the second 
trapezoid* which Is evaluated and summed Into FLDZ. Each successive 
grid bar Integral closes another trapezoid for the transverse Integral 
such that the final grid bar closes the last trapezoid and FLDZ 
contains the result of the aperture Integration for the theta polarl- 
zatlon of the field. FLDZ Is stored In FIELDZ(I) and the computation 
continues with the next set of angles. When the Integrals have been 
evaluated for all angles In the pattern, the subroutine returns to 
the main procedure. 

3.7 Main Procedure 

The main procedure orchestrates the principal subroutines 
described earlier In the section to make a radiation pattern calcula- 
tlon. The variables Introduced are defined below: 

ET0T(2,400) - 

A summing array for the complex values returned by INTGRT. 
FIELDY and FIELDZ are summed into ET0T(1,*) and ET0T(2,*)* 
respectively. 

YFLD(75) , ZFLD( 75 ) - 

Arrays containing the magnitudes of the unnormal Ized $- 
and G-polarlzed electric field strengths, respectively. 

In a final antenna pattern. 
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I 

I 

j PWER ( 75 ) - 

An array containing the unnormalized power for angles In 
a final antenna pattern. 

§ 

S PR(2,500) - 

|* An array containing the perimeter points of all panels of the 

| reflector for use by the APRMAP subroutine. The y- and 

z- coordinates are denoted by the first subscript being one or 
two, respectively. Panels are punctuated in the array by a 
vary large number In the y-coordlnate location. 

* 

MtVL - 

The stack pointer for the PR array. It denotes the next 
location pair [PR(l,-)» PR(2,-)3 to receive a perimeter point 

I YLO, YHI - 

The lowest and highest y-coordlnate possessed by any panel 
yet processed. 

ZLO, ZHI - 

The lowest and highest z-coordlnate possessed by any panel 
yet processed. 

I SUM - 

The Index of the last pair of locations In the ETOT array to 
receive or supply pattern data. 

NANG - 

The number of angles In the present pattern. 

FMAXY, PMAXZ - 

The maximum and o-polarlzed electric field magnitudes, 
respectively. In a particular pattern. 

PWR - 

The total power corresponding to FMAXY and FMAXZ. 

0 - 

m 

1 

I 


The current value of the Incremented angle In a particular 
pattern. 
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FMYDB, FMZOB, PWRMDB - 

The decibel values of FMAXY, FMAXZ, and PWR, respectively. 

DBZZ, DBZY - 

The relative power In dB of the e-polarlzed electric field 
with respect to FMZOB and FMYDB, respectively. 

OBYZ, OBYY - 

The relative power In dB of the ^-polarized electric field 
with respect to FMZOB and FMYDB, respecively. 

The procedure begins by defining the decibel function and calling 
NPUT to read the reflector antenna parameters and computational 
requests. A loop Is then begun which Is performed once for each panel. 
APRTUR Is called to create the Illumination array and trace the rays to 
the aperture plane; then QUANTZ Is called to prepare the aperture for 
Integration. APRPLT Is called If the output control code dictates. 

The dimensions of the aperture are checked with the current maximum and 
minimum values (YLO, YHI, ZLO, and ZHI), after which the perimeter 
points are stored In the PR array In locations Indicated by the stack 
pointer MLVL. MLVL Is updated and punctuation for that panel Is stored 
after Its perimeter points. The punctuation Is used by APRMAP to dis- 
tinguish between separate panels. 

Patterns are stacked In the ETOT array as shown In Figure 13. 
INTGRT Is called once for each pattern request. Each time a pattern Is 
returned by INTGRT, FIELOY and FIELOZ are summed Into the portion of 
the stack containing that pattern. ISUM Is the pointer Indicating the 
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last pair of locations in ETOT to receive values. It Is updated after 
every pattern and reset to zero after every panel. 

After these steps have been performed, APRTUR is called to begin 
computation of the next panel. The process is repeated until all 
panels are finished, after which an APRMAP plot is generated, if 
specified by the output option code. 

The rest of the main procedure is concerned with printing tables 
of electric field strengths for each pattern request. The first 
pattern is taken from the ETOT array and the magnitudes of the and 
e- polarized fields are pieced in the YFLD and ZFLD arrays, respec- 
tively. The maximum and minimum of each of these arrays is found and 
assigned to FMAXY and FMAXZ, respectively. If these values are 
greater than 10"^, they are converted to dB and the pattern from 
which they are taken is computed, tabulated, and plotted. If one 
or both of them are less than 10"^, the entire pattern is set to 
-60 dB. This is an indication to the programmer that the magnitude 
of that polarization of the field is insignificantly small or that 
some problem exists in the calculation. Once the pattern maximum has 

been deemed large enough to output the pattern, the lowest value 

-15 

converted to dB for an individual angle is 10 . Any value smaller 

than this is assigned the value of -100 dB. This provides 100 dB 
of dynamic range for the weakest pattern qualifying for computation. 
The above precautions are exercised to prevent problems with the 
logarithm function in the comouter. 

The relative field values (DBZZ, DBZY, DRYZ, and DBYY) are 
computed and the table is printed one line at a time. PL0T4 is 


called to make a printer plot of the table entries. The above steps 
are repeated until all requested patterns have been printed. 
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4. OUTPUT CONTROL AND PROGRAM OPTIONS 
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4.1 Introduction 

In order to better monitor the calculation, a large quantity of 
output Is available from the program. Selection of specific output Is 
accomplished by means of a three-digit output option code which Is 
supplied to the program along with the antenna parameters In the Input 
file. This section describes the printing, plotting, and permanent 
storage abilities of the program and how the code may be used to exer- 
cise them. 

The output option specification requires between one and eight 
cards, depending on which values are supplied for the output option 
code. These cards appear just before the pattern request cards in the 
input file. The option code, known to the computer as N0PT(3), is 
punched in the first three columns of the first card. Another field on 
the first card is for the variable NLIST. The remaining seven cards, 
if required, are for the dimensioned variable IIIST(IOO). These 
variables will be explained later in the section. 


4.2 Printing Options 

An example of the printed output from REFLECTR may be found in 
Appendix E. The portion of the output prior to the illumination data 
is not under control of the option code. The panel illumination and 
quantizing data (also in Appendix E) are under the control of the first 
digit of the code. If it is desired to print panel data for all panels, 
N0PT(1) is specified as two. It may be desirable to print this infor- 
mation only for several specific panels. In this case, N0PT(1) is 
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specified as one and the total number of panels for which data Is 
requested Is entered for the NLIST parameter. NLIST Is right-justified 
In column 10 of the first card. If N0PT(1) Is one, the program will 
look for the numbers of NLIST specific panels on the following cards 
according to a 1615 format. Panels are numbered by the program In the 
order In which they are received. If It were required to print panel 
data for panels 6, 12, and 19, for example, NLIST would be given the 
value 3 and the first three 15 fields on the next card would contain 
the numbers 6, 12, and 19. These need not be supplied In any specific 
order and no more option cards would be required. If no panel data Is 
to be printed, N0PT(1) is specified as zero. Printed output following 
the panel quantizing data; that is, the computed radiation pattern of 
the antenna, is not under control of the code. 

4.3 Plotting Options 

Two types of plotted output are available from REFLECTR: a plot 
of data points in the aperture plane (APRPLT) and a plot of the panel 
aperture boundaries on the aperture plane (APRMAP). APRPLT plots the 
aperture data points for a single panel (an example appears in Appendix 
E). APRMAP, which also appears in Appendix E, plots the boundaries of 
all of the apertures in one plot. Plotting is controlled by the second 
digit of the output option code, N0PT(2). If all of the panel apertures 
are to be plotted, N0PT(2) is specified as two. If only specific panel 
apertures are to be plotted, N0PT(2) is specified as one. In this case, 
the procedure followed is identical to that for the printing option. 

If both the printing and plotting options are specified as one, the 
panels to be printed and plotted must be the same. That is, only one 


set of cards for ILIST follow the option code card; thus, there Is no 
provision for printing and plotting different sets of panels. APRMAP 
Is called automatically If any panels are plotted. If no plotting is 
to occur, N0PT(2) Is specified as zero. Note that NLIST and ILIST only 
need to be provided If N0PT(1) or N0PT(2) equals one. 

4.4 Disk Storage Options 

REFLECTR has the ability to write the post-QUANTZ P array to a 
disk file for later examination or to be retrieved at the end of the 
program. This option Is controlled by the third digit of the code, 
N0PT(3). If N0PT(3) Is specified as two; the P array will be stored for 
all panels, then retrieved and printed at the end of the program. If 
N0PT(3) Is one, the P array will be stored for all panels. Finally, 
specifying N0PT(3) as zero will cause the program not to store any 
aperture data. The data may be examined at a later time by submitting 
a job which calls the PTLIST subroutine described In Section 5.4. 


5. THE SUPPORT CODE 


5.1 Introduction 

The support subroutines are subroutines which do not participate 
directly In the radiation pattern calculation. Instead, they aid the 
calculation In some way--e1ther by preparing output or by performing 
simple data processing tasks. These subroutines do not appear In 
Appendix C with the principal subroutines. Instead, a listing of each 
Is provided along with Its description in this section. 

The plotting package does not appear in this section. Plotting 
subroutines tend to be extremely installation-dependent and It was felt 
that a plotting package using local software would not be useful In 
general application. Agrawal's version of APRPLT [3] makes a good 
example on which to base a new version. There Is enough Information in 
Section 3.7 and the test case plots In Appendix E with which to 
synthesize a new version of APRMAP. 

5.2 PL0T4 

PL0T4 is the printer plot subroutine which produces the output 
shown in Figure 14. It plots the table of relative field strengths 
comprising ths radiation pattern over a 60 dB dynamic range with one 
angular increment to a line and one dB to a print column. If a value 
to be plotted is more than 60 dB below the maximum value, a "<" charac- 
ter is printed at -60 dB. The angular axis is labelled every five 
angular increments and the lower plot border is always extended to the 
next even multiple of five increments. A listing of the subroutine 
appears in Table 1. 
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TABLE I. PL0T4 


SUBROUTINE PLOT4 (NAME, FMAX, P, NT, MINOR, AMINOR) 

IMPLICIT REAL*8 (A-H,0-2) 

REALMS NAME (8) , MINOR 
DIMENSION LINE (71) 

DIMENSION F (NT) , AMINOR (3) 

DATA IBLANK, LESS, I DOT, IDASH, I/1H , 1H<, 1H . , 1H-, 1HI/ 

PRINT 10, NAME, MINOR 

10 FORMAT (///25X,8A8,///7X,A5, * (DEC) -60' ,7X, '-50',7X, '-40',7X, 

•-30‘,7X,'-20',7X,'-10’,7X, , 0',9X,'+10',/20X, , *',9X,'*',»X, 
' ,9X, ' , 9X , ' * ' , 9X , • * * , 9X , ' * ' , 9X , ' * ' ) 

DO 20 J-1,70 

20 LINE (J) -IDASH 

DO 22 K-l ,71,5 
22 LINE (K ) «I 

DO 170 J ■ i , NT 
IP (J.EQ.l) GO TO 41 
DO 29 N-1,70 

29 LINE (N) bIBLANK 

TESTbDFLOAT(J-1)/5.0 

IP ( (TEST-DINT (TEST) ) .EQ. 0.0) LINE (2)-IDASH 
IP (LINE (2) .EQ. IDASH) LINE (70)-IDASH 
41 PLACE-F (J ) -PMAX+60. 0 

IP (PLACE) 90,80,70 
70 LINE(1)bI 

K-PLACE+1.5 
LINE(K)bIDOT 
00 TO 99 

80 LINE (1 ) “IDOT 

00 TO 98 

90 LINE (1) -LESS 

98 LINE (2 ) “IBLANK 

99 PRINT 105, LINE 

10S FORMAT (20X , 7 1A1 ) 

ANGLE-AMINOR (1)+ (J-l)*AMINOR (3) 

IF (LINE (70) . EQ. IDASH) PRINT 149, ANGLE 
149 FORMAT ('+' ,8X,F8. 3) 

170 CONTINUE 

DO 515 L-2,70 

515 LINE (L)-IBLANK 

LINE (1)-I 
II-NT+1 
.J.I-NT+fi 

00 520 L-II.JJ 
TRST»DFL0AT (L-l ) /5. 0 

t F ( (TEST-DINT (TEST) ) . EQ. 0. 0) GO TO 524 
520 PIUNT 105, LINE 

524 no 525 M * 1 , 70 

525 LINE (M) «IDASH 

DO 545 N-1,71,5 

545 LINE (N ) » I 

PRINT 105, LINE 
ANGLE»ANGLE4-5 . * AM I NOR ( 1' 

PRINT 149, ANGLE 

RETURN 

END 
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The lines of the plot are composed one at a time in the character 
array LINE(71). The array begins and ends with the character "I". On 
every fifth multiple of the angular Increment, LINE(2) and LINE(70) are 
loaded with the character and the value of the angle Is printed 
beside the line. The dot Is placed by using the value of the relative 
field strength to compute the Index of LINE which will locate It 
properly with respect to the dB axis. The upper and lower plot borders 
are generated by loading all 71 locations with dashes, then writing an 
"I" over every fifth dash. 

The calling arguments are NAME(-), FMAX, F(-)» NT, MINOR, and 
AMIN0R( • ) . NAME(*) Is a l-by-8 array of eight-byte variables containing 
the title of the plot. The title must be 64 characters long or padded 
with blanks to make 64 characters. FMAX Is the maximum unnormalized 
field strength magnitude for the entire pattern In dB. F( • ) Is an 
array containing the unnormalized magnitudes of the pattern field 
strengths In dB. NT Is the number of angles in the pattern. Finally, 
MINOR Is the name of the angle to be varied in the pattern, and 
AMIN0R(1) is for the low, high, and incremental values of this angle 
for 1 equal to one, two, and three, respectively. 

5.3 10 

10 is the output control function subroutine. A referral to the 
output control variables is made by supplying 10 with an integer code 
for the particular output operation and the number of the panel cur- 
rently under computation. 10 returns a 1 if the operation Is permitted 
and a 0 if it is not. 
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A listing of 10 Is presented In Table 2. INTENT Is the output 
operation code variable. The codes are: 

1. Print panel Illumination and quantizing data. 

2. Plot the P array after quantizing (APRPLT). 

3. Create a map of panel projections with the plotter (APRMAP). 

4. Write the P array to a disk file after quantizing. 

5. Read the P array from the disk file and print It at the end 
of the program. 

ITER, or Iteration, Is the number of the present panel. 

5.4 Utility Routines 

Utility routines are those subroutines which perform standard data 
processing tasks. Included In this section are SETM, MOVEM, INTPL, 
MULT32, PTSORT, and PTLIST. All of these routines are relatively 
simple; thus, only a brief description of each will be provided. 

SETM Is listed In Table 3. It loads N locations of a singly- 
dimensioned array (A) with the constant X. Because of the way arrays 
are Indexed In the computer, SETM may also be used to load X Into rows 
of doubly-dimensioned arrays--such as the P array. For example, the 
Nth row of the P array could be loaded with X using the statement 
CALL SETM (X,P(1 ,N) ,5) . 

SETM Is used primarily In this capacity. 

MOVEM appears In Table 4. Its function Is to copy N values from 
one singly-dimensioned array (AOLD) Into another array of the same 
type (ANEW). Like SETM, MOVEM may also be used to operate on rows of 
doubly-dimensioned arrays. The main use of MOVEM in this program Is to 
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TABLE 2. 10 


FUNCTION XO (INTENT, ITCH) 

IMPLICIT RUL*« (A-H,0-t) 

INTEGER 9 UK FCC 

COMMON /CONTRL /NOFT (3) ,NLXST, XOFT, XLXST (100) 
00 TO (20,31,40,50,60) ,INTtMT 
20 IF (NOFTtl) ,EQ.O) CO TO tO 

XF (NOFT (1 ) ,EQ« 2) GO TO 91 

22 DO 25 I«1 ,NLI9T 

IF ( I LIST (X ) • tO* ITER) GO TO 91 
25 CONTINUE 

GO TO 90 


30 

IF 

(NOFT (2) 

.err. 

0) 

CO 

TO 

91 


CO 

TO 90 






40 

IF 

(NOFT (2) 

.to* 

0) 

GO 

TO 

90 


IF 

(NOFT (2) 

.to. 

J) 

CO 

TO 

91 


CO 

TO 22 






50 

IF 

(NOPT (3) 

*CE. 

11 

GO 

TO 

91 


CO 

TO 90 
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IF 

(NOPT (3) 

• EQ. 

2) 

CO 

TO 

91 


90 10*0 

RETURN 


91 I0«1 

RETURN 

END 


TABLE 3. SETM 


SUBROUTINE SETH (X, ANEW, N) 
IMPLICIT REAL’S (A~H,0-2) 
DIMENSION ANEW (N) 

IP (N . LE. 0) RETURN 
DO 100 1*1, N 
ANEW ( I )*X 
100 CONTINUE 
RETURN 
END 


TABLE 4. MOVEM 


SUBROUTINE MOVEM (AO LD, ANEW, N) 
IMPLICIT REAL’S (A-H,0-2) 
DIMENSION AOLD (N ) , ANEW (N ) 

I F (N , LE« 0) RETURN 
DO 100 1*1, N 
ANEW ( I) *AOLD (I ) 

100 CONTINUE 
RETURN 
END 
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copy l-by-5 vectors such as PBLK(S) Into rows of the P array, as shown 
below: 

CALL MOVEM (PBLK, P(1,N),5). 

INTPL is a linear interpolation subroutine used for the edge point 
Interpolations In the aperture plane. It Is listed In Table 5. PL0(*) 

and PHI ( - ) are arrays containing the data for the two known points with 
the lower and higher z-coordlnates , respectively. Z(*) Is an array 
containing the z-coordlnates of the endpoints for the grid bar on which 
the Interpolation Is to be performed. The z-coordinate of the endpoint 
between the known points Is assigned to the Interpolated point PINT(*)< 
The phase and both components of the electric field are Interpolated 
based on the relationship between the z-coordlnate of the endpoint and 
those of the known points. 

MULT32 premultiplies a 3-by-2 matrix (OPB) with a 3-by-3 matrix 
(OPA) and stores the result In a 3-by-2 array (ROP). It Is listed In 
Table 6. 

PTSORT sorts a 5-by-n array A(5,n) In ascending order, first with 
respect to the numbers stored In the first column A(1,n), then with 
respect to the second column, A(2,n). The first and second columns 
correspond to the y- and z-coordinates of the data points In the 
P array, respectively. A listing of PTSORT is presented In Table 7. 
NROW is the length of the array in the second dimension. SWAP(*) Is an 
Intermediate array used to store information during a row exchange. 

The operation of the subroutine is simple. The first row Is compared 
to every other row in the array. If any row is found which should 
appear before the first row, the information in the rows Is 
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TABLE 5. INTPL 


SUtftOtrilfE INTPL (PLO, PN I , PINT, 2) 
IMPLICIT REAL*! (A-H # 0 -2) 

DIMENSION PLO (5) ,PHI (5) , PINT (5) ,2(2) 
20 * 2 ( 1 ) 

IP (PLO (2) • GT. 2(1)) 2C*2 (2 ) 

PAC- (2C-PLO (2 ) ) / (PHI (2) -PLO (2) ) 

DO 10 1-1,5 

10 PINT (I ) *PLO (I ) ♦FAC* (PHI (I) -PLO (I)) 

RETURN 
END 


TABLE 6. MULT 32 


SUBROUTINE MULT 3 2 (ROP,OPA,QPB) 

IMPLICIT REAL*# (A-H,0-2) 

DIMENSION ROP (3,2) ,OPA(3, 3) , CP# (3, 2) 

RQP ( 1 , 2) -OPA ( 1 , 1) *OPB ( 1 , 2) fOPA ( 1 , 2) *OPB ( 2, 2) +OPA (1,1) *OP# (3,2) 
ROP (2,2) *OPA ( 2 , l) *OPS ( 1 , 2) +OPA (2,2) *OP# (2, 2)*0PA (2, 3) *OP# (3,2) 
ROP (3,2) *OPA (3,1) *OPB ( 1 , 2) *OPA (1,2) *OP# (2, 2) ♦OPA (3,3) *OP# (3 , 2) 
ROP (1 , 1 ) *OPA ( 1 , 1 ) *OPB (1,1) ♦OPA (1,2) *0P1 (2,1) ♦OPA (1,3) *OPB (3,1) 
ROP (2,1) *0PA (2 , 1 ) *OPB (1,1) ♦OPA (2,2) *OP# (2,1) +OPA (2,3) *OPB (3, 1 ) 
ROP(3,l)*OPA(3,l)*OPB(l, 1)+0PA(3,2)*0PB(2,1)40PA(3,1)*0PB(3 # 1) 
RETURN 
END 


TABLE 7 PTSORT 


SUBROUTINE PTSORT (A, NCOL, NROW) 
IMPLICIT REAL*# (A-H ,0-2 ) 

DIMENSION A (NCOL, NROW) , SWAP (5 ) 

DO 900 1-2, NROW 
DO 900 .1*1, NROW 

IF { A(1,J)-A(!,I-1) ) 100,200,900 
100 DO ISO K » I , NCOL 

SWAP (K ) • A ( K , I - 1 ) 

A(K, 1-1 )«A(K,.l) 

ISO A (K , J ) *SWAP (K ) 

CO TO 900 

200 IF ( A(2,J)-A(2, I-i) ) 220,900,900 
220 DO 250 K*1 , NCOL 

SWAP (K)*A(K,I-1) 

A (K, 1-1 ) *A (K , J) 

250 A (K , J ) *SWAP (K ) 

900 CONTINUE 
RETURN 
END 
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Interchanged. After scanning the entire array, the first row will be 
that which has the lowest numbers In the first and second columns. The 
second row Is then compared to the third and greater rows, and so forth, 
until the last row Is reached. At this point, the array Is sorted. 

PTLIST reads the aperture data points which have been stored on 
disk earlier In the program and prints them In tables, one for each 
panel. NLIST Is the number of points in the aperture and IPNL is the 
number of the panel corresponding to the aperture. A negative number 
read in the NLIST field indicates the end of the file. A listing of 
PTLIST appears in Table 8. 
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TABLE 8. PTLIST 


SUBROUTINE PTLIST 
KEAL*8 0(5) 

READ (7,18) NLIST, IPNL 

IF (NLIST. LE.O) GO TO 8 

PRINT 20, NLIST, IPNL 

DO 4 1-1, NLIST 

READ (7, 21) Q 

PRINT 22, 0 

CONTINUE 

GO TO 2 

PRINT 19 

STOP 

FORMAT (218) 

FORMAT (///10X, '********* END OF RECORD ON DATA SET 7 ***•****••) 
FORMAT ('l',//l OX, 'LISTING OF APERTURE PLANE POINTS AS SUPPLIED TO 
.RADPATV/10X, 16, * POINTS INTE®ATED FOR PANEL NUMBER', 14, 

.///14X, 'Y -COORDINATE ' ,8X, ' Z -COORDINATE ', 10X, 'EY (V/M) ' ,12X, 

. 'EZ (V/M) ',1 OX, 'PHASE (RADIANS)'//) 

FORMAT (5D16. 10) 

FORMAT (10X,F15.8,6X,F1S.8,SX,D1S. 10,4X,D1S. 10, 4X, FIS. 8) 

RETURN 

END 
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6. EXAMPLE AND FIRST TEST CASE 
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6.1 Introduction 

The use of the program Mill now be demonstrated. In addition* a 
card-by-card description of the general Input file will be provided 
before considering the specific antenna configurations. The example 
Is taken from work reported by Rusch [5] and was used In an earlier 
report on REFLECTR prior to the work reported herein. The reflector 
will be described In general, then analyzed by the program in three 
roughly equivalent configurations. In Case I, the single-panel option 
will be used to represent the reflector as a single piece of material 
with an elliptical edge. In Case II, the reflector will be represented 
as a single panel with the edge defined by 20 perimeter points. This 
will closely approximate the elliptical edge. Finally, the reflector 
will be described as three panels such that these panels together form 
the surface of Case III. The object of the procedure is to compute 
nearly identical patterns in all three cases and to demonstrate agree- 
ment with the calculations in [5]. The first case verifies that the 
point processing changes made in QIIANTZ will not adversely affect the 
performance of the program. The second case shows that a reflector 
edge may be adequately represented by a polygon. The third representa- 
tion shows that it is numerically justified to superpose the patterns 
of the component panels to obtain the total pattern and that the 
algorithm is working correctly. In each case, the complete input file 
will be developed and listed. 


I ' ii „• i.1 i -t! i*! 1 1 ! . ' r !! : ' ■N!ii»«iiil ! ! ; !f n ! TH t; Hf '1hKMSH(!»l|jUt!fJitf! 


The antenna, show) in Figure 15, is an offset-fed parabolic dish 
with a focal length of 31.5 inches. Measured in the foces-centered 
reference coordinates , the physical reflector lies between 21.8 and 
84.1 degrees in the theta coordinate. The aperture projects onto the 
yz plane as a one-meter circle. The phase center of the feed antenna 
is located at 

x = -1.854 in. 
y = 0.0 
z = -3.966 in. 

in the reference system with the horn tilted upward 40 degrees. 

The feed system origin is placed at the phase center of the feed 
antenna. The system is inclined 40 degrees in order to make 6 * 90°, 

$ * 180° correspond to the bores ight of the antenna. This coordinate 
system rotation is given by ALPHA = 0°, BETA = 0°, and GAWiA = -40°. 
The translation vector must be expressed in feed system coordinates. 
This is found by 

x 1 = 1.854 sin 40° - 3.966 sin 40° = -1.129 in. 

y' = 0.0 

z' = 3.966 cos 40° + 1.854 sin 40° = 4.230 in. 

The aperture plane should be placed just forward of the reflector 
surface. Some clearance should be allowed to insure that the over- 
illumination rays encounter the aperture plane before striking the 
reflector surface. The maximum x-coordinate of the surface in this 
example may be found by using the polar form of the parabola equation. 


.-JL. 



FIG. 15. OFFSET- FED PARABOLA 
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- . 2f 
1 7 cos 6, 


max 


2(31.5) 

1 + cos 66. 


yr » 


45.94 In. 


x ■ -r cos 0_. w ■ -45.94 cos 68.2° ■ -17.06 In. 

Ifl&X 


The aperture plane will be placed at XC ■ -16.0 In. to provide the 
necessary clearance. 

The feed antenna for this reflector Is a scalar horn with a 15 dB 
amplitude taper at 34.1 degrees. Rusch [5] computed field strengths 
for this antenna In one-degree Increments off bores Ight In both the E- 
and H-planes. These values are used by the FILL routine described In 
Section 3.4 to assign field strengths to the Illumination array. The 
radiation Is ♦'-polarized with a wavelength of 2.367 inches [4.99 GHz], 

6.3 General Input File 

To help clarify the organization of the input file for REFLECTR, 
the general Input file is explained below. The three examples In the 
following sections will provide files for actual computations. 


Card(s) Information 

1-4 Title cards 


5 FEEDO-3), ALPHA, BETA, GAMMA, XLAM 

6 SURFCE, NPNL, AORORF, BELLP, PSI 

7 PLNPNTO-3), PLN0RM(l-3) 

8 XC, YCBL, ZCBL, HFMABL, HFMIBL 


9-N Any data required by the FILL subroutine. 

(N+1 )-(N+4) Four extreme points (X, Y, Z) on the edge of the 
reflector. One point goes on each card (single 
panel option only). 
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Card(s) Information 

N+5 NOPT, NLIST 

(If NOPT specifies that only certain panels are to 
be printed or plotted, cards containing the list of 
these panels follow this card.) 

N+6 MAJOR. AMAJOR, MINOR, AMIN0R(1-3) 

(Pattern request cards. There will be between one 
and five of these cards In this location— the order 
of the requests is not Important.) 

N+7 DONE typed In the first four columns of the card. 


These cards are followed by the panel data unless the single panel 
option Is exercised. The organization of the panel data appears below. 


Card(s) 

1 

2- ( NPERIM+1 ) 
NPERIM+2 


Information 

NPERIM, SURFCE, NPTPPL 

(X, Y, Z) Perimeter points. 

(One perimeter point Is put on each card.) 

AORORF or AORORF, 8ELLP or AORORF, PSI or PLNPNT(l-3), 
PLN0RM(l-3) depending on which parameters are needed 
to describe the surface specified by SURFCE. 


The files for the individual panels are concatenated In the main input 
file after the DONE card. The order in which the panels appear Is 
unimportant; however, they will be numbered by the program In this 
order. 


6.4 First Configuration 

Case I uses the single-panel option. In addition to the general 
information computed in the previous section, four points defining the 
aperture boundary ellipse must be found. Since these points are used 
to find the maximum and minimum 9' and <j»' values for the illumination 
array, these four points must represent the extreme points on the dish 
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In each angular coordinate. The 0' extrema are the uppermost and lower- 
most points on the dish and may be found as follows: 

Upper point - 


r . 2f . 2(31,5) 

1 + cos B 1 + cos 68.2° 

max 

x ■ -r cos e mmw * -17.060 in. 
max 

y » 0.0 

z » r sin * 42.654 In. 
max 


45.94 In. 


Lower point - 


r 

x 


2f 


i + cos 6 


min 


r-r^srl^ =31.58 in. 


-r cos 0 . = -31.416 in. 


y = 0.0 

z * r sin 6^ s 3.246 in. 


The two points representing the y-extrema of the dish are almost 
exactly the extrema as well. The z-coordlnates of these points are 
Identical . 


z 


min 


+ Z max 
1 


42.654 + 3.246 
2 


22.95 in. 


The dish is one meter wide and symmetrical with respect to the xz plane 


y 


. 39.37 ... 

- * 


±19.685 in. 
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Finally, the paraboloid equation provides the x-coordlnate. 
x ■ ^ -f ■ -24.244 In. 

Table 9 shows the Input file for Case I. The first four cards 
contain title Information which Is reproduced at the beginning of the 
printout. The coordinate system Information (FEED, ALPHA, BETA, and 
GAMMA) and the wavelength (XLAM) appear on Card 5. Cards 6 and 7 con- 
tain Information about the surface of the reflector. Card 6 Is for 
SURFCE, NPNL, AORORF, BELLP, and PSI. Card 7 Is for PLNPNT and PLNORM. 
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For this reflector, SURFCE, NPNL, and AORORF are given values of 4, 0, 
and 31.5, respectively. None of the other parameters are required for 
this surface, so all are given the value zero. Card 8 Is for aperture 
plane Information— XC, YCBL, ZCBL, HFMABL, and HFMIBL. In this case, 
there Is no aperture blockage, so HFMABL and HFMIBL are both set to 
zero. Cards 9 through 46 are the tables of electric field strengths 
for the particular feed antenna In one-degree Increments. The values 
begin at 0 degrees and progress to 90 degrees as one reads text— line 
by line from left to right. The 41= 180® plane follows the e = 90° plana 
The card order that follows is peculiar to the single-panel option. 

Cards 47 through 50 are for the extreme points on the dish. Unlike the 
perimeter points, they may be entered In any order. Card 51 Is for the 
output option code. Here the computer Is Instructed to exercise the 
print and plot options for all panels and to write the aperture plane 
point array onto a disk file at the end of the QUANTZ subroutine. Since 
It is not requested that output options be exercised for any particular 
panels, NLIST has been left blank and there are no panel number cards 
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I 

r 

r 

* for HIST following Card 51. Cards 52 and 53 art tht radiation pattarn 

|" requests. Patterns are requested In the * ■ 0° plane for e from 68.0* 

to 98.0° by 0.5°, and In the 6 ■ 83* plane for + from -15.0* to 15.0* 
j by 0.5*. The next card has DONE typed In the first four columns, which 

signifies that no more pattern requests follow. All of the above cards 
*• are read by the NPUT subroutine. This completes the Input file for 

]' Case I. The output appears in Appendix E. 
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6.5 Second Configuration 

In Case II, the program will run normally with 20 points on the 
dish edge selected as perimeter points. It was previously noted that 
the dish edge projects to a one-meter circle on the yz plane. The 
perimeter points are calculated by choosing twenty equlspaced points on 
this circle and solving the paraboloid equation for the appropriate 
x-coordinate. The Input file for this test case Is shown In Table 10. 
Th^ only changes In the file before Card 47 are the title cards and the 
NPNL. parameter. NPNL Is changed from 0 to 1 to Indicate that the panel 
Is to be treated normally and perimeter points are to be expected. The 
four extreme points used In Case I are deleted so that the option code 
card becomes Card 47 and the pattern requests become Cards 48 and 49. 
Card 50 Is the DONE card. The panel data begins on Card 51 with NPERIM, 
SURFCE, and NPTPPL. NPTPPL Is entered as 2250 to match the default 
value used in the single-panel option. The next NPERIM cards each con- 
tain the coordinates of a perimeter point (x, y, z) In Its order of 
appearance as one traces around the dish edge. The final card In the 
panel data file is reserved for wahtever parameters are required for the 
panel surface. In this case, the focal length of the paraboloid is 
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TABLE 10. CASE I INPUT 


OfFStT-fO PARABOLIC DISH 
XLLUK2JUTX0N XNTtBFOLATtD 
PABfL ONt* 20 PUTS PAJtfl 


mavumm 

TAIL**; *0 BLOCKAOt) 


PASSU - 

TWO- 21 NTI PAWSL TRItlt- 21 PWTS 


CASS Tims 


1 

1 

HAY 22, 

1960 NCSU 

POHAi 60T0U 

pans 

HXLLS80N0 

I 

5 

-1.129 

0.0 *.229 

0.0 

0.0 

-60.0 

2.367 


A 

A 

3 31.5 

0.0 

0.0 




7 

0.0 

0.0 0.0 

0.0 

0.0 

0.0 


r 

6 

• 16.0 

0.0 0.0 

0.0 

0.0 



I 

9 

1.000000 

0.99651* 

0.99*251 


0.967122 

0.977230 

* 

10 

0.96*6*5 

0.9*9*70 

0.931627 


0.911655 

0.669713 


11 

0.665567 

0.639603 

0.612006 


0.762963 

0.752735 

r - 

12 

0.721*7* 

0.669*09 

0.656752 


0.623711 

0.590*66 

| 

13 

0.557261 

0.52*260 

0.491659 


0.459590 

0.426226 

I. 

1* 

0.597712 

0.366164 

0.3397*7 


0.312513 

0.266573 


IS 

0.262009 

0.236690 

0.217275 


0.197210 

0.176732 


16 

0.16106* 

0.1*6621 

0.133000 


0.120967 

0.1105*6 

| 

17 

0.101622 

0.09*129 

0.067955 


0.062956 

0.078977 

1 

11 

0.075636 

0.073357 

0.07137* 


0.069737 

0.066321 

* 

19 

0.067025 

0.065772 

0.06*510 


0.063205 

0.061636 


20 

0.06039* 

0.056675 

0.057260 


0.055609 

0.053666 

? - 

21 

0.052066 

0.050215 

0.0*6333 


0.0*6*36 

0.04*5*7 

i 

i 

22 

0.0*2671 

0.0*0616 

0.036990 


0.037167 

0.035*15 

i - 

23 

0.033663 

0.03200* 

0.030391 


0.02665* 

0.027391 


26 

0.025996 

0.02*666 

0.023366 


0.02216* 

0.021002 


25 

0.019913 

0.016906 

0.017961 


0.017126 

0.016331 

} 

26 

0.015576 

0.01*656 

0.01*177 


0.013550 

0.012969 

1 

27 

0.012*96 






* ■■ 

26 

1.000000 

0.996*66 

0.99*139 


0.966670 

0.97679* 


29 

0 963975 

0.9*6527 

0.930566 


0.91029* 

0.66762* 

f * 

30 

0.6633*6 

0.637063 

0.809177 


0.779691 

0.7*9*26 

1 

31 

0.716007 

0.6656*9 

0.653170 


0.620162 

0.567093 

i . 

32 

0.55*100 

0.521397 

0.469156 


0.*575** 

0.426712 


33 

0.396795 

0.367920 

0.3*018* 


0.313662 

0.266*90 


36 

0.26*666 

0.2*2263 

0.221307 


0.201620 

0.183606 

j 

35 

0.167258 

0.152156 

0.136*73 


0.126161 

0.115165 


36 

0.105*17 

0.096637 

0.08933* 


0.062606 

0.0771*9 


37 

0.0722*8 

0.067993 

0.06*276 


0.061002 

0.056077 


38 

0.055*27 

0.052967 

0.050705 


0.0405*3 

0.046*70 

? *■ 

39 

0.0*4*65 

0.0*2513 

0.0*0607 


0.0387*2 

0.036916 

1 

60 

0.035129 

0.033363 

0.031662 


0.030026 

0.026*26 

4 * 

61 

0.026678 

0.025369 

0.023961 


0.02259* 

0.021291 


62 

0.020050 

0.01667* 

0.017760 


0.016709 

0.015719 


63 

0.01*790 

0.013916 

0.01310? 


0.012336 

0.011623 

i 

66 

0.01095* 

0.010329 

0.009**6 


0.009203 

0.006696 


65 

0.006223 

0.007760 

0.007366 


0.006976 

0.006609 


66 

0.006265 







67 

221 







66 

PHI 

0.0 THtTA 

66.0 

98.0 

0.5 


§ 

69 

THETA 

63.O PHI 

-15.0 

15.0 

0.5 


v* 

50 

DONE 
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TABLE 10. (continued). 


47 221 

1(8 PHI 0.0 THETA 68.0 98.0 0.5 

49 THETA 83.0 PHI -15.0 15.0 0.5 

50 DONE 


51 

20 

4 2250 


52 

-24.244 

19.685 

22.950 

53 

-22.028 

18.722 

29.033 

54 

-20.029 

15.925 

34.521 

55 

-18.443 

11.571 

38.875 

56 

-17.424 

6.083 

41.672 

57 

-17.073 

0.0 

42.635 

58 

-17.424 

-6,083 

41.672 

59 

•18.443 

-11.571 

38.875 

60 

-20.029 

-15.925 

34.521 

61 

-22.028 

-18.722 

29.033 

62 

-24.244 

-19.685 

22.950 

63 

-26.460 

-18.722 

16.867 

64 

-28.459 

-15.925 

11.379 

65 

-30.046 

-11.571 

7.025 

66 

-31.064 

-6.083 

4.228 

67 

-31.415 

0.0 

3.265 

68 

-31.064 

6.083 

4.228 

69 

-30.046 

11.571 

7.025 

70 

-28.459 

15.925 

11.379 

71 

-26.460 

18.722 

16.867 

72 

31.5 




needed (AORORF). Were the panel pert of an ellipsoid, both AORORF and 
BELLP would be expected on this card. A parameter received In a panel 
data file replaces any value supplied previously. The output received 
from this configuration Is shown In Appendix E. 

6.6 Third Configuration 

Case III requires the synthesis of three separate panels from the 
single parabolic dish. The division of the surface and the location of 
the perimeter points are shown in Figure 16. The perimeter points on 

the edge of the dish are the same as those used in Case II. The points 

on the face of the dish were found by relating the y- and z-coordinates 
with line equations and using the paraboloid equation to find the 
corresponding x-coordi nates. The density of illumination Is kept close 

to that of the first two cases by requesting that each panel be Illumi- 

nated by 800 rays. 

The input file is shown in Table 11. It is identical to the file 
for Case II through Card 50 except for the title cards, the NPNL speci- 
fications, and the NPTPPL parameter. The illustrative purposes, the 
panels have been input in the order shown in Figure 16, although the 
panels may be input in any desired sequence. Organization of the panel 
input files follows from Case II. The output appears in Appendix E. 

6.7 Results 

The results of these calculations are shown in Figures 17(a,b). 

The solid line represents Rusch's calculations. The circles, triangles, 
and squares represent the results of Cases I, II, and III, respectively. 
Whenever the Case II or III points fall on those computed from Case I, 
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TABLE II. CASE m INPUT 


1 OFFSET-FED PARABOLIC DISH RUSCH/ AEROJET 

2 ILLUMINATION INTERPOLATED FROM TABLES; NO BLOCKAOE 

3 REFLECTOR IS ONE PANEL DEFINED BY 20 PERIMETER POINTS 


CASE TOO 


4 

HAT 22, 

1980 NCSU 

PGMR: BOTULA 

FCLTT : 

HILLSBORO 

5 

-1.129 

0.0 4.229 

0.0 

0.0 

-40.0 

2.367 

6 

A 

1 31.5 

0.0 

0.0 



7 

0.0 

0.0 0.0 

0.0 

0.0 

0.0 


6 

-16.0 

0.0 0.0 

0.0 

0.0 



9 

1.000000 

0.998514 

0.994251 


0.987122 

0.977230 

10 

0.964645 

0.949470 

0.931827 


0.911855 

0.889713 

11 

0.865567 

0.8:9603 

0.612006 


0.782983 

0.752735 

12 

0.721474 

0.689s 09 

0.656752 


0.623711 

0.590488 

13 

0.557281 

0.524260 

0.491659 


0.459590 

0.428226 

14 

0.397712 

0.366184 

0.339747 


0.312513 

0.286573 

15 

0.262009 

0.238890 

0.217275 


0.197210 

0.178732 

16 

0.161864 

0.146621 

0. 133000 


0.120987 

0.110546 

17 

0.101622 

0.094129 

0.087955 


0.082958 

0.078977 

18 

0.075836 

0.073357 

0.071374 


0.069737 

0.068321 

19 

0.067025 

0.065772 

0.064510 


0.063205 

0.061836 

20 

0.060394 

0.058875 

0.057280 


0.055609 

0.053868 

21 

0.052066 

0.050215 

0.048333 


0.046438 

0.044547 

22 

0.042671 

0.040818 

0.036990 


0.037187 

0.035415 

23 

0.033683 

0.032004 

0.030391 


0.028854 

0.027391 

24 

0.025998 

0.024666 

0.023388 


0.022164 

0.021002 

25 

0.019913 

0.018906 

0.017981 


0.017128 

0.016331 

26 

0.015576 

0.014856 

0.014177 


0.013550 

0.012989 

27 

0.012496 






28 

1.000000 

0.998486 

0.994139 


0.986870 

0.976794 

29 

0.963975 

0.948527 

0.930586 


0.910294 

0.887824 

30 

0.863348 

0.837063 

0.809177 


0.779891 

0.749428 

31 

0.718007 

0.665649 

0.653170 


0.620182 

0.587093 

32 

0.554100 

0.521397 

0.489158 


0.457544 

0.426712 

33 

0.396795 

0.367920 

0.340184 


0.313682 

0.288490 

34 

0.264668 

0.242263 

0.221307 


0.201820 

0.183806 

35 

0.167258 

0.152158 

0.138473 


0.126161 

0.115165 

36 

0.105417 

0.096837 

0.089334 


0.082808 

0.077149 

37 

0.072248 

0.067993 

0.064278 


0.061002 

0.058077 

38 

0.055427 

0.052987 

0.050705 


0.048543 

0.046470 

39 

0.044465 

0.042513 

0.040607 


0.038742 

0.036916 

40 

0.035129 

0.033383 

0.031682 


0.030028 

0.028426 

41 

0.026878 

0.025389 

0.023961 


0.022594 

0.021291 

42 

0.020050 

0.018874 

0.017760 


0.016709 

0.015719 

43 

0.014790 

0.013918 

0.013102 


0.012338 

0.011623 

44 

0.010954 

0.010329 

0.009746 


0.009203 

0.008696 

4S 

0.008223 

0.007780 

0.007366 


0.006976 

0.006609 

46 

0.006265 
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TABLE 

II. 

(continued) 


51 

20 

4 800 


I 

52 

- 27.320 

0.0 

22.950 

I, 

53 

- 26.428 

- 1.653 

25.225 


54 

- 25.411 

- 3.306 

27.500 

1 

55 

- 24.269 

- 4.959 

29.775 

L 

56 

- 23.001 

- 6.612 

32.050 


57 

- 21.607 

- 8.265 

34.325 


58 

- 20.088 

- 9.918 

36.600 

1. 

59 

- 18.443 

- 11.571 

38.875 

60 

- 20.029 

- 15.925 

34.521 


61 

- 22.028 

- 18.722 

29.033 

i 

62 

- 24.244 

- 19.685 

22.950 

1 

63 

- 26.460 

- 18.722 

16.867 


64 

- 28.459 

- 15.925 

11.379 

r 

i. 

65 

- 30.046 

- 11.571 

7.025 

66 

- 30.033 

- 9.918 

9.300 

67 

- 29.895 

- 8.265 

11.575 


68 

- 29.631 

- 6.612 

13.850 

j 

69 

- 29.241 

- 4.959 

16.125 

1 

* 

70 

- 28.726 

- 3.306 

18.400 


71 

- 28.086 

- 1.653 

20.675 

i 

72 

31.5 



i 

i 

73 

21 

4 800 



74 

- 27.320 

0.0 

22.950 

1 - 

75 

- 27.257 

2.812 

22.950 

i 

i 

76 

- 27.064 

5.624 

22.950 

77 

- 26.755 

8.436 

22.950 

t - 

78 

- 26.316 

11.249 

22.950 

§ 

i 

79 

- 25.751 

14.061 

22.950 


80 

- 25.060 

16.873 

22.950 


81 

- 24.244 

19.685 

22.950 

i 

82 

- 22.028 

18.722 

29.033 

i 

m - 

83 

- 20.029 

15.925 

34.521 


84 

- 18.443 

11.571 

38.875 

r 

85 

- 17.424 

6.083 

41.672 

86 

- 17.073 

0.0 

42.635 

A - 

87 

- 17.424 

- 6.083 

41.672 


88 

- 18.443 

- 11.571 

38.875 

=9 * 
1 

89 

- 20.088 

- 9.918 

36.600 

•# - 

90 

- 21.607 

- 8.265 

34.325 


91 

- 23.001 

- 6.612 

32.050 

r 

92 

- 24.269 

- 4.959 

29.775 

1 

93 

- 25.411 

- 3.306 

27.500 


94 

- 26.428 

- 1.653 

25.225 

T 

95 

31.5 




I 


jta 





ET 












83 


TABLE II. (continued). 


96 

21 

4 800 


97 

- 27.320 

0.0 

22.650 

98 

- 28.086 

- 1.653 

20.675 

99 

- 28.726 

- 3.306 

18.400 

100 

- 29.241 

- 4.959 

16.125 

101 

- 29.631 

- 6.612 

13.850 

102 

- 29.895 

- 8.265 

11.575 

103 

- 30.033 

- 9.918 

9.300 

104 

- 30.046 

- 11.571 

7.025 

105 

- 31.064 

- 6.083 

4.228 

106 

- 31.415 

0.0 

3.265 

107 

- 31.064 

6.083 

4.228 

108 

- 30.046 

11.571 

7.025 

109 

- 28.459 

15.925 

11.379 

110 

- 26.460 

18.722 

16.867 

111 

- 24.244 

19.685 

22.950 

112 

- 25.060 

16.873 

22.950 

113 

- 25.751 

14.061 

22.950 

114 

- 26.316 

11.249 

22.950 

115 

- 26.755 

8.436 

22.950 

116 

- 27.064 

5.624 

22.950 

117 

- 27.257 

2.812 

22.950 

118 

31.5 




mm 
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only the circles are plotted. In all cases, the H-plane results agree 
almost perfectly with Rusch's calculations. The major deviation In the 
E-plane pattern occurs where Rusch predicts a beam shoulder and 
REFLECTR predicts a shallow null. This Is probably due to the fact 
that such points on the beam are extremely dependent on the placement 
of the feed phase center. Rusch accounted for the difference In phase 
center between the E- and H-plane feed patterns, while we did not. 
Still, the E-plane patterns agree well with Rusch's computed pattern. 
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7. SEGMENTED SPHERICAL REFLECTOR TEST CASE 

7.1 Description of the Antenna 

In order to test the algorithm with measured data, a segmented 
spherical reflector antenna was constructed and tested by NASA at the 
Langley Research Center. The reflector surface consisted of 54 
triangular-shaped planar facets, the vertices of which were located on 
a sphere with a 24-Inch radius. The center of the surface was elevated 
21 degrees above the negative x-axis. The boundary of the reflector 
was hexagonal in shape and measured 11.6 inches from point to point. 

The surface was fabricated from flberglas molded onto a machined metal 
negative and sprayed with conducting paint. The feed antenna was a 
dual mode horn with a 16 dB amplitude taper at 23 degrees. Illumi- 
nating radiation was y-polarized at a frequency of 35 GHz. The center 
of the horn aperture was located at -12.52 inches on the negative 
x-axis and the feed antenna was elevated 40 degrees above horizontal. 

A diagram of the antenna appears In Figure 18. 

7.2 Computation and Results 

For computational purposes, the phase center of the feed horn was 
assumed to be 0.5 (a) or 0.169 inches behind the aperture. The aperture 
plane was placed at x ■ -18.7 inches. The triangular panels were 
represented by three perimeter points corresponding to their three 
vertex points as computed for the mold surface. Each panel was Illumi- 
nated with 150 rays. Patterns were computed in the e = 91° plane with 
4> varied from 0 to 16 degrees and In the $ * 0° plane with e varied 

\ 
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from 81 to 101 degrees. The Incremental angle was 0.5 degrees In both 
cuts. The algorithm performed the computation In 6.3 minutes. 

The results of the computation are compared with the measured 
patterns In Figures 19(a,b). The H-plane pattern compares extremely 
well with the measured results. The E-plane pattern exhibits a slightly 
wider beamwldth and a higher amplitude of the fourth sldelobe. The 
fourth sldelobe Is a grating lobe caused by the periodicity of the 
reflector surface In the y-coordlnate. It Is present regardless of 
the polarization of the Illuminating radiation. It should be noted 
that the program exactly predicts the location of the lobe as well as 
the fact that It Is more steeply sloped on Its Inward edge. The level 
of the lobe, however. Is about 6 dB lower than predicted. 

One possible explanation for these discrepancies Is that the 
vertex points used In the program were computed for the mold and not 
measured from the actual surface. The reflector may have sprung out 
somewhat in the xy plane after being removed from the mold, which would 
cause the reflector surface analyzed by the program to be different 
from that used in the pattern measurements. The exact location of the 
vertex points on the flberglas reflector will be measured In the future 
and the computation will be repeated in hopes of obtaining better 
agreement In the E-plane. 





FIG. 19b. SPHERE E-PLANE PATTERN 
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8. CONCLUSIONS 

An algorithm has been developed to compute radiation patterns for 
large segmented reflector antennas by extending an existing algorithm 
for smooth-edged, single-panel reflectors. Included In this algorithm 
are new techniques for aperture boundary definition, feed antenna des- 
cription, output formatting, ano reflector Illumination. Equations for 
parabolic cylinders have been Included so that the program can be used 
to analyze deployable umbrella reflector antennas. 

The algorithm has been tested by comparing calculations with 
calculations made by Rusch [5], and with measured data, with the results 
presented In Sections 6 and 7. 

Because of the discrepancies between the calculated and measured 
E-plane patterns of the faceted spherical reflector. It Is recommended 
that the reflector surface be measured and the program rerun to see If 
the pattern might be better predicted. If the spherical test case 
does not achieve better results, a return to the theory underlying the 
program would be required W determine what electromagnetic effects are 
contributing to the t ror and how they might be included in the program 
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APPENDIX A 

DEVELOPMENT OF THE APERTURE PLANE INTEGRALS 


In a source-free region, the electric field at a point Is given 


by 

( 1 ) 

where sinusoidal time variation is assumed. 

Supplying the Lorentz gauge condition, this can be rewritten as 


| = Y1 • j 

jwey 


Ju), 


Ivx 


F. 


( 2 ) 


By the equivalence principle [ 6 ], the tangential aperture electric 
field becomes a magnetic current source of strength M = -2(n x E). The 
radiated electric field strength is thus 


where 


-*■ 1 

E = V x F 




(3) 

(4) 


The integration over the source volume is independent of the 
observation coordinates acted on oy the v operator. ( 3 ) and ( 4 ) are 
combined and the curl is taken of the integrand. 


. +-jkr 

E = -47 /' v x (Ml—) dV* . 

To expand the integrand, we make use of the vector identity 
v x = (v$) xNH(’xfi) 

and note that M is a function only of source coordinates, thus 


(5) 

( 6 ) 
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4>( V x M) « 0 


-jkr 


V* = 7 (~ ) 

■ ♦ } Utn + j itiilj 


e -jkr , 

f (jk + 1). 


If ft is represented as 


ft - xM x + yMy + zM z , 


the cross product is 
-jkr 


vxftl 


fr.M 


(z-z*r 

A 

X 

L 2 

r y 

r 


+ l-M 


(x-x')1 

A 

y > 

1 X 

L A 

r z 

r 


+ l-M 

l*-*!)- ♦ M 

(y-y-Il 

A 

z 

l - y 

r x 

r 



- ikr 


M ■ 


We now make the far field approximations: 


(A) 

(B) 


[»♦*]* H 

r * r Q - [x* sinecos* + y ’ sinesin^* + z'cose 
(See Figure Al) 


(C) x-x' s x 

y-y' = y 

z-z' = z 

(D) -■ = sinecos^ 
^ = sines in<j> 
7 = cose 


( 7 ) 

(8) 

( 9 ) 


Using these approximations in the integrand ( 9 ), we obtain 
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Ofe- jkr - 


|^-M z s1nes1n<|> + M y cose ft 

j^-M x cose + MjSlnecos^ y 

£-M y sinecos<l> + M x sinesin4> z 


exp{ j k[x * sinecos 4 > + y'sinesin<|>+ z'cose]}. 


The electric field integral becomes 

t m -jk<j' jkr 0 j Gdv , 

4irr 0 v 1 


where 


G = -M sines in<|>+ M cose 
* y 


-M x cose + M z sinecoS(« y 


r-v 


necos<t»+ M x sinesin4j z 


• exp {jk[x'sinecos(j> + y'sinesin<|> + z'cose]}. 

The source is an aperture in the y'z' plane, ft is given by 
ft « -2(n x t) 

- 2[-x x (E x ft + E^ + E 2 z)] 

■ -2z A E y + 2 pE z . 


M x = 0 

M y = 2E z 
H 2 = -2E y 

and the electric field in the far field is giver, by 
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. J [E «1ne*1n* + E,ca 5 ae Jk k’s , "« , "*« , «se] d > dz 

x 4ltr 0 s' y 2 


. -J2>e' j>r ° i [.£ „s1necosde jk[y ' s1n9s '"* +2 ' cos6 Vdz' 


4irr 


s' 


(14) 

(15) 


. -j 2 ke - Jkr 0 j [. E s11)9cosJ ,4kCy's1ne5ln^'co88] dy , dl , - 

^ * q s 1 ^ 

A A 

The vectors <l> and e in spherical coordinates are given by 

A A A 

<f> * -x sini}> + y cos<j> 

A A A A 

e * x cos^cose + y sin^cose - z sine. 

-2kp"j kr n 

E 0 = — jhnr / ( U exp{ j k [y ' s i n e s i n 4> + z 'cose]}dy 'dz ' 


(16) 

(17) 

(18) 

(19) 


where 


U = cos<Kose[EySinesin(j> + E z cose] + sin^cose [-E y sinecos4»3 
+ sine[E z sinecos4>] 

= E (sinesin^cos^cose - sinesin 4 >cos<(>cose) 

J 

2 2 

+ E_(cos ecos4> + sin ecos«(>) 


U = E z cos<j> 


S -Jkr f 


E^ = - —4^ - -/V exp { j k [y ' s i ne s i n<j> + z'cose]} dy'dz' 


( 20 ) 

( 21 ) 


where 


V = -sin4>[E sinosin<j> + E z cose] + cos$[-E s i necos<j >3 

J J 

2 2 

- E y (-sin (j>sine - cos <j>sine) - E z s i n^cose 


V = -E sine - E z sin4>coso . 

Jr 


(22) 
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We are Interested In the relative magnitudes and phases of these 
field components. Furthermore, the relative phase of each component 
needs only to be preserved with respect to Itself. Hence, the constant 
premultlpller may be dropped. The Initial phase of the fields on the 
aperture Is accounted for as * In the phase term. We have 

E 0 * // E 2 cos<p exp{j[ky's1nes1n<t> + kz'cose -*]}dy'dz' (23) 

Aperture 

E, ■ //[E w s1ne + E_coses1n$ Jexp{j[ky'sines1n<j> + kz'cose -*]}dy'dz' 

* y Z (24) 

Aperture 

These integrals are evaluated numerically in the program. 
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APPENDIX B 

DEVELOPMENT OF RAY BACKTRACKING EQUATIONS 

The object of tracing a ray back to the feed antenna Is to express 
a point on the reflector surface by its angles e' and <t>' In the feed 
spherical coordinate system. 

Let us define a new coordinate system which Is parallel to the 
feed system but centered on the reference origin. These coordinates 
are denoted x", y" , and 2 ''. The point on the reflector Is expressed 
as a vector and rotated into the double-prime system by the Inverse of 
the original rotation matrix. 




X 


X' ' 

A' 1 


y 

5 

y' ' 



z 


_ 2 " _ 


The translation to the feed system is accomplished by adding the 
translation vector O'O. The x, y, and z components of this vector are 
called FEED(l), FEED ( 2 ) , and FEED(3) , respectively. 

The distance R of the point from the feed origin is found by 
taking the magnitude of the vector 0M3 + V" = V'. 

R ■ t(x ') 2 ♦ (y'l 2 ♦ (z ') 2 ] ,/2 
The angle e' is found by 
z‘ = R cos e ' 
e* = cos" 1 ^-) 


The angle <p' is found by 




y' * R sin 6' sin <f>' 

♦' ■ s1n '' ( ft sin e' > • 
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R and sin e 1 will always be positive. If y 1 Is positive, the desired 
<j>' Is less than it radians. The arcsine function will return an angle 
between 0 and ir/2. If y' Is zero, the desired ' Is ir, but the arcsine 
returns 0. Finally, If y' Is negative, the desired <j>' Is greater than 
ir, and the arcsine returns a negative angle. It Is concluded that If 
the returned angle Is subtracted from it, the desired angle Is obtained. 

♦' ' * - s1n " 1 *R""sfrr r * 
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APPENDIX C 

ADDITION OF PARABOLIC CYLINDERS 

Standard form of the parabolic cylinder equation shall be with the 
directrix in the xz plane, opening In the positive x direction. The 
focal point of the directrix lies on the coordinate origin. This sur- 
face is shown In Figure Cl. 



FI6. Cl. STANDARD POSITION 

The equation for the surface in this position is 

z 2 = 4f(f + x) (1) 

where f is the focal length of the directrix. We now allow the surface 
to be rotated about the x-axis. Rotation is considered as positive 
when counter-clockwise if the observer looks In the negatlve-x direction. 
The angle of rotation is denoted by i>. The coordinate system rotation 
is pictured in Figure C2. 
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FIG. C2. COORDINATE ROTATION 

An arbitrary point P may be rotated from the unprimed to the primed 
coordinates by the relations 

z' * z cos <j< - y sin (2) 

y' * z sin \p + y cos v . (3) 

The parametric equations for a ray are 
x - RB^ - B ]? 

y * RB 21 - b ?2 (4) 

z » RB 31 - B 32 

First account for the rotation of the surface by substitution of z' In 
(2) for z In (1). 

(z') 2 « 4f(f + x) 

(z cos ^ - y sin ij») ■ 4f(f + x) 
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2 2 2 2 
z cos i|> - 2 yz sin * cos * + y sin i|> » 4f(f ♦ x). (5) 

Then substitute the parametric expressions for x, y, and z. 

2 2 

(RB 31 - R 32 ) cos ij/ - 2 (RB 21 - B 22 )(RB 31 - B 32 )cos 0 ln(j» 

♦ (RB 21 - B 22 ) 2 s1n 2 4> ■ 4f(f ♦ RB^ - B 12 ) 

(R 2 b|, - «B 31 B 32 ♦ B^)cos 2 * - 2(R 2 B 21 B 31 -R^ 1 B 22 -RB 21 B 32 +B 22 B 32 ) 
cos^slmj* + (R 2 b | 3 - 2RB 21 B 32 + B^jsln 2 * - 4f 2 -4fRB 11 

+ 4fB 12 « 0. (6) 


The last equation (6) Is of the form 
AR 2 + BR + C * 0 


where 


2 2 2 2 
A - B 31 cos i|> - 2B 21 B 31 cos^s in\|< 4 B 21 s1n * 


p 

B * ~2B 3 -|B 32 COS 4* ^ 2 i B 22 4 B 2 ^ 6 32 )co$i^s1n<p 

2 

- 2B 21 By 2 s1n <|< - AfBjj 


(7) 

( 8 ) 

(9) 


2 2 2 2 2 
C « B 32 cos > - 2B 22 B 32 cos ^ s1n '*' + B| 2 s1n> - 4f c 4 4fB^ 2 (10) 

Expressions (8), (9), and (10) are evaluated by the program and 
(7) Is solved to find the Intersection point of the ray with the 
surface. 

Using (5), the inside normal of "he surface may be found at any 
point by taking the negative of the gradient. 

-v$(x,y,z) « n(x ,y,z) (11) 


where 


2 2 2 2 
$(x,y,z) ■ * cos - 2yz sln^cos^ ♦ y sin * - 4f(f 4 x) ■ 0 
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* 4f 

2 

■ 2z sln^cos* - Zy sin 

2 

-v$ 2 ■ 2y slftycos* - Zz cos tf> . 


02 ) 

03) 

04) 


A unit vector Is obtained by normalizing this vector. The common 
factor of 2 Is removed from each component prior to normalization such 
that 

n ■ x 2f/NMAG ♦ y(z sln^cos^ - y sin^|<)/W1AG 
+ z(y slm^cosiji - z cos >)/NMAG 

where 

NMAG ■ [4f^ + (z s in^cos^ - y s1n 2 <j>) 2 + (y sln^cos^ - z cos 2 *) 2 }^ 2 
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PROGRAM REFLBCTR 

IMPLICIT REAL*8 (A-H.O-Z) 

REAL*8 MAJOR (5) .MINOR (5) 

INTEGER SURPCE 

COMPLEX*16 ETOT (2,400), FIELDY (400) ,PIBLDZ(400) 
COMMON/APRPRM/NPTPPL,NPERIM 
COMMON/CONTRL/NOPT (3) ,NLIST, IOPT, ILIST (1 00) 

COMMON/D IMENS/YDIM , ZDIM , YCT, ZCT 
COMMON/EXTENT/YM IN , YMAX , ZM IN , ZMAX 
COMMON/MATH/P I , PI 2 , PID2 , DTOR , RTOD 

COMMON/PATTRN/ETOT,AMINOR (3,5) ,AMAJOR(5) , MINOR, MAJOR, NANGLE (5) 
COMMON/PARAMS/AORORP, BELLP, PSI , PLNPNT (3) ,PLNORM(3) ,PEED(3) , 

. ALPHA , BETA , GAMMA , XLAM , XC , SURPCE , NPNL ,NP0INT 

DIMENSION P (5,2750) , YPLD (75) ,ZPLD (75) ,PWER (75) ,PR (2,500) 

DATA DONE/5HDONE /,MLVL,NPARTS/0, 7/ 

DATA YLO, YHI , ZLO, ZHI/1 . 0D+10, -1 . 0D+10, 1 . 0D+10, -1 . 004-10/ 

FDB (X) *20. 0*DLOG10 (X) 

PDB(X)«10.0*DLOG10(X) 

MAXPTS-2750 
CALL NPUT (P , NPAT) 

DO 400 1*1, NPNL 

CALL APRTUR (P , I ) 

PRINT 777 

CALL QUANTZ (P,NPERIM,I) 

PRINT 778 

IF (10(3,1). EQ. 0) GO TO 80 
ISW*I 

IF (IOPT. EQ* 1 ) ISW— 1 
CALL APRPLT (P ,NPOINT, ISW) 

PRINT 780 
80 CONTINUE 

IF (YMIN.LT.YLO) YLO-YMIN 
IF (YMAX. GT. YHI) YHI*YMAX 
IF (ZMIN.LT.ZLO) ZLO-ZMIN 
IF (ZMAX. CT.ZHI) ZHI *ZMAX 
DO 95 L=i ,NPERIM 
PR ( l , L+MLVL) *P (1,L+N POINT) 

95 PR (2 , L+MLVL) “P (?. , L+N POINT) 

mlvl*mlvl+nper:m-h 
PR ( 1 , MLVL) -1 . 0D+40 
IS UM*0 

DO 200 K*1,NPAT 

CALL INTGRT (P, MAJOR (K) ,AMAJOR (K) , AMINOR (1 , K) , FIELDY, PIELDZ) 
PRINT 779 
NANfUNANGLF, (K) 

DO 150 t. = 1 , NANG 

ETOT (1 ,L+TSUM) -ETOT (1 ,L+ISUM) +FIELDY (L) 

150 ETOT (2, L+ISUM) =ETOT (2 , L+ISUM) +FIELDZ (L) 

200 I S UM - 1 S UM +NAN G 

400 CONTINUE 
PRINT 781 

IF (IOPT.EQ.l) GO TO 420 

IF (10(2,1). EQ. 0) GO TO 420 

YDIM-YHI-YLO 

Z DIM*ZHI -Z LO 

YCT = (YHI+YLO) /2.0 

ZCT- (ZHI + ZL0)/2.0 

CALL APRMAP (PR, NPNL, -1) 

PRINT 782 

420 ISUM-U Vl^CrPi-- * -j • .. 

DO 770 1*1, NPAT np v/ ' , * ' 

NAN G*N ANGLE (I) " • 
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FMAXY— 1.0D+40 

FMAXZ— 1.0D+40 

DO 450 J>1,NANG 

YFLD (J ) "CDABS (ETOT (1 , J+XSUH) ) 

ZFLD (J ) “CDABS (ETOT (2, J+ISUM) ) 

FMAXY-DMAX1 (FMAXY, YFLD (J) ) 

450 FMAXZ >DMAX1 (FMAXZ , ZFLD (J ) ) 

ISUM-ISUM+NANG 
D»AMINOR (1,1) 

FMYDB- -60. ODO 

FMZDB- -60. 0D0 

PWRM DB— 60.000 

PWR-FMAXZ+FMAXZ+PMAXY+PMAXY 

IP (FMAXY. CT. 1.00-10) FMYDB-FDB (FMAXY) 

IP (FMAXZ. GT. 1.0D-10) PMZDB-PDB (FMAXZ) 

IP ( PWR .CT. 1.00-10) PWRMDB>PDB( PWR ) 

PRINT 600, MAJOR (l ) , AMAJOR (I) .MINOR (I) , (AMINOR (J, I) ,J-1,3) 

600 FORMAT (1H1 , ///24X, 

. 'TABLE OF ELECTRIC FIELD STRENGTHS (DB) ',/'+', 23X, 

' 

. //7T§XT' PRINCIPAL PLANE OF CUT IS ' ,Ji57 r - ',P8.3,' DEG' 

//19X, 'ANGLE '.AS,' FROM',P8.3,' TO',P8.3,‘ BY',P6.3,' DEG') 
PRINT 666, MINOR (I) 

666 FORMAT (//13X,A5,4X,'DB(Z/Z) ' ,4X, ' DB (Y/Z) • ,4X, • DB (Z/Y) * , 5X, 

. ' DB ( Y/Y ) ' , 5X, ' PWRDB ' ,/) 

DO 700 K-l.NANG 

PKER (K ) -PWRMDB -1 00 . 0D0 

DBY -FMYDB -100. 0D0 

DBZ «PMZDB -100. ODO 

PWR-ZFLD (K ) *ZFLD (K ) +YFLD (K ) *YFLD (K ) 

IF (YPLD (K ). GT. 1 . OD-15) DBY-FDB (YFLD (Kl ) 

IF (ZFLD(K) .GT.l. OD-15) DBZ-FDB (ZFLD (K ) > 

IF (PWR.GT. 1.00-30) PWER (K) -PDB (PWR) 

IF (FMYDB . EQ. -60 . ODO) DBY— 60. ODO 

IF (FMZDB. EO. -60. ODO) DBZ— 60. ODO 

DB7.Z-DBZ-FMZDB 

DBYY-DBY-FMYDB 

D3ZY -DBZ -FMYDB 

DBYZ-DBY-FMZD3 

PWRDB “PWER (K ) -PWRMDB 

PRINT 690, D, 03?. Z , DRYZ , DBZY , DBYY , PWRDB 
690 FORMAT (10X,F9. 3,5F11.5) 

D-D+AMINOR (3,1) 

YFLD (K ) -DPY 
ZFLD (K ) -DBZ 
700 CONTINUE 

PRINT 750, FMAXZ, FMZDB, FMAXY, FMYDB 
750 FORMAT (//15X, 'MAXIMUM FIELD VALUES-' //I 5X, 

. ’ 20 LOG (MAX (FIELD-Z))-20LOG(’,1PE15.7,')-', 0PFI2 , 7//1 5X , 

20 LOG (MAX (FIELD-Y))“20LOG( , ,1PE15.7,')>',OPP12.7) 

PRINT 755, NPARTS 
755 FORMAT (// 1 4X , 

. ' INTERPOLATION NIWBFR USED FOR INTEGRATION IS ',15) 

PRINT 76 5 , MAJOR (I ) , AMAJOR (I ) 

765 FORMAT(1H1,///20X, 'PRINCIPAL PLANE - ’,A5,P7.3,' DEGREES') 

CALL PL0T4 (64H NORMALIZED Z-COMPONENT OF SECONDARY PATTERN (DB) 

, FMZDB, ZFLD, NANG, MINOR (I) , AMINOR (1,1)) 

PRINT 765, MAJOR (1) , AMAJOR (I) 

CALL PL0T4 (64H NORMALIZED Y-COMPONENT OP SECONDARY PATTERN (DB) 

, FMYDB, YFLD, NANG, MINOR (I ) , AMINOR (1,1)) 

PRINT 765, MAJOR (I ) .AMAJOR (I) 

CALL PLOT4 (64H NORMALIZED POWER PATTERN (DB) 

. , PWRMDB, PWER, NANG, MINOR (I) ,AMINOR(l,I) ) 
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770 


775 

776 

777 

778 

779 

780 

781 

782 


CONTINUE 

IF (10(4,1) .80. 1) WRITE (7 ,775) 
IF (10(5,1). EC- 0) STOP 


REWIND 7 
CALL PTL1ST 
FORMAT ( ' -1234') 

FORMAT (/' FINISHED INPUT 

FORMAT (/' FINISHED APRTUR 

FORMAT (/' FINISHED "'UANTZ 

FORMAT ( ’ FINISHED INTGRT 

FORMAT ( • ***** EXECUTED APRPLT ***** 

FORMAT (/• PATTERN COMPUTATIONS COMPLETE 

FORMAT (/' ***** EXECUTED APRMAP ***** 

STOP 

END 


’) 

') 

*) 

') 

') 

') 


SUBROUTINE NPUT(P.NPAT) 

IMPLICIT REAL*8 (A-H.O-Z) 

REAL'S MAJOR (5) .MINOR (5) 

C0MPLEX*16 ETOT (2 ,400 ) 

INTEGER SURFCE 

COMMON/BLOCKG/YCBL, ZCBL, HFMABL, HFMIBL 
COMMON/CONTRL/NOPT (3) ,NLIST, IOPT, ILIST(IOO) 

COMMON/FEED/EP (91) , ET (91) ,NP,NT 
COMMON /MATH/P I , PI 2 , PID2 , DTOR , RTOD 

COMMON/PARAMS/AORORF , BE LLP, PS I , PLNPNT (3 > ,PLN0RM(3) ,FEED{3) , 

ALPHA , BETA , GAMMA , XLAM , XC , S URFCE, NPNL , NPOINT 
COMMON/PATTRN/ETOT, AMINOR (3, 5) .AMAJOR (5 ) , MINOR , MAJOR, NANGLE (5) 
DIMENSION P (5, 2750), TITLE (40) 

DATA DONE/5HDONE / 

IOPT-O 

READ 200, TITLE 

READ (1 , 10) FEED, ALPHA, BETA, GAMMA, XLAM 
10 FORMAT (7F10. 4) 

READ ( 1 , 20 ) .5 URFCE , NPNL , AORORF , RELLP, PS I , PLNPNT, PLNORM 
20 FORMAT (II, 14X, 15, 3F10. 4/6F10. 4 ) 

READ (1 , 30) XC.YCBL, ZCBL, HFMABL, HFMIBL 
30 FORMAT (5F 10 . 4 ) 

READ ( 1 , 40) EP 
READ (4,40) ET 
40 FORMAT (5F15. 6) 

IF (NPNL. LE . 0) READ il , 50) ( (P (I , J) , I *1 , 3) , J-l , 4 ) 

50 FORMAT (3F1 0.4) 

READ (1,60) NOPT.NLIST 
60 FORMAT (31 1 , 2X , I 5) 

IF (NOPT(l) .EQ, l.OR.NOPT (2) . EQ. 1) READ(1,70) (I LIST (I ) , 1-1 ,NLIST) 
70 FORMAT (1615) 

ISUM-0 

NPAT-1 

77 HEAD (1,80) MAJOR (NPAT) .AMAJOR (NPAT) .MINOR (NPAT) , (AMINOR (I,NPAT) , 

. 1-1,3) 

80 FORMAT(A5,5X,F10.4,A5,5X, 3F10.4) 

IF (MAJOR (NPAT) . EQ. DONE) GO TO 98 

NANGLE (NPAT)- (AMINOR (2 , NPAT) -AMINOR (1 , NPAT) ) /AMINOR (3, NPAT) *1. 5 
IF (MANGLE (NPAT) .GT. 75) GO TO 88 
IS IW-ISUM+N ANGLE (NPAT) 

NPAT-NPAT+1 

IF (NPAT. LT. 6) GO TO 77 

PRINT 573 

STOP 


i 


no 


83 

98 

105 

108 

120 

130 

140 

150 

160 

170 

180 


190 

200 

573 

574 

575 

576 

577 

578 


579 


PRINT 574 
STOP 

IP <ISUH,LE,400) GO TO 105 

PRINT 575# I S UM 

STOP 

NPAT»NPAT-i 

DO 108 L«1,ISUM 

ETOT (1 , t) - (0 . 0D0, 0. 0D0) 

ETOT (2 , L) ■ (0 . 0D0 1 0. 000) 

GOTO (120,130, 140,150,160) ,SURFCE 

PRINT 578, TITLE, XLAM, FEED, ALPHA, BETA , GfcMMA, PLNPNT, PLNORM 
GO TO 170 

PRINT 579, TITLE, XLAM, PEED, ALPHA, BETA, GMtMA, AORORF, BE LLP 
GO TO 170 

PRINT 580, TITLE, XLAM , FEED, ALPHA , BETA , GAMMA , AORORF 
GO TO 170 

PRINT 581, TITLE, XLAM, FEED, ALPHA, BETA, GAMMA, AORORF 
GO TO 170 

PRINT 582, TITLE, XLAM, FEED, ALPHA, BETA, GAMMA, AORORF, PS I 

IF (NPNL.GE.l) GO TO 180 

IOPT-1 

NPNL*1 

PRINT 583, ( (P (I , J) , I»1 , 3) , J-l , 4) 

PRINT 584, XC , YCBL, ZCBL , HFMABL , HFMIBL,NPNL 

PRINT 586 

PRINT 587 

PRINT 590, EP 

PRINT 586 

PRINT 588 


PRINT 590, ET 
PRINT 576.NPAT 
DO 190 M«1,NPAT 

PRINT 577, MAJOR (M) , AMAJOR (M) , MINOR (M) , (AMINOR (KK, M) , KK«1 ,3) 
FORMAT (10A8) 

FORMAT (' ********** ERROR - MORE THAN 5 PATTERN ', 
•CALCULATIONS REQUESTED **•*****•*•) 

FORMAT ( ' ********** ERROR - MORE THAN 75 ANGLES IN', 

• ONE PATTERN REQUEST ***********) 

FORMAT ( * ********** ERROR - REQUESTED' , 1 5, ' ANGLES TO BE ', 


•CALCULATED EXCEEDS AVAIL. STORAGE **********•) 

FORMAT(// 

' NUMBER OF PATTERN GROUPS REQUESTED ’,15/) 


FORMAT (5X , A5 , • - • ,F10. 4 , 10X, A5, ' FROM • ,F10. 4, • TO*,F10.4,' 8Y', 


F10.4) 

FORMATflHl,////, 13X, 'PLANAR REFLECTOR FAR FIELD RADIATION V 

• PATTERN CALCULATION*///' ',10A8/' ',10A8 
/' ' , 10AB/' ' , 10A8// 

• INPUT PARAMETERS- ' ,// 

• WAVELENGTH OF ELECTRIC FIELD ',F9.4/ 

' LOCATION OF COORDINATE ORIGIN WRT FEED (X,Y,Z) ',3F8.3 

/' FEED ROTATION ANGLES (ALPHA, BETA, GAMMA) ',3F8.3 

/' A POINT ON THE REFLECTOR SURFACE (X,Y,Z) '.3F8.3 

/' COMPONENTS OF UNIT NORMAL TO SUFACE (X,Y,Z). ',3F8.3 

/) 


FORMAT (1H1,////,11X, 'ELLIPTICAL REFLECTOR FAR FIELD RADIATION •/ 
• PATTERN CALCULATION'///' ',10A8/' *,10A8 


/' ' , 10A8/' ' , 10A8// 

' INPUT PARAMETERS- ' ,// 

' WAVELENGTH OF ELECTRIC FIELD ',F9.4/ 

' LOCATION OF COORDINATE ORIGIN WRT FEED (X,Y,Z) ',3F8.3 

/' FEED ROTATION ANGLES (ALPHA, BETA, GAMMA) ',3r8.3 

/' MAJOR AXIS OF THE ELLIPTICAL REFLECTOR ',P78.3 


/' MINOR AXIS OF THE ELLIPTICAL REFLECTOR ',f8.3) 


m 


580 


581 


582 


583 


584 


586 

587 

588 
590 


FORMAT (1H1 #//// ,11X, 'SPHERICAL REFLECTOR FAR FIELD RADIATION */ 

/• mom/- mom// '*"“"“““»’>««•///• mom/, mom 

' INPUT PARAMETERS- , .. 

' WAVELENGTH OP ELECTRIC FIELD tpa a/ 

' LOCATION OF COORDINATE ORIGIN WRT FEED (X. Y, 2) ! ! " " • * 1P8 3 

FEED ROTATION ANGLES (ALPHA, BETA, GAMMA) .' 3M *3 

RADIUS OF REFLECTOR SPHERE . 3 ? 

FORMAT (1H1,////,UX, ‘PARABOLIC REFLECTOR FAR FIE LD *RADIATION ' *•/ 

/• MOM/' MOM// MOCOUtlO..///. MOM/' MOM 

* INPUT PARAMETERS- , /y 

#// 


/' 

/* 


/* 


WAVELENGTH OF ELECTRIC FIELD •'»« a, 

LOCATION OF COORDINATE ORIGIN WRT FEED (X , Y,*Z j ! .’ ! ! ! • ! 3F8 3 
FEED ROTATION ANGLES - ' ' #3F8.3 


/’ FOCAL LENGTH OF THE REFLECTOR^* GMWA) •'* 0 * 3 ? 

FORMAT (1H1,////,HX, *,ARMOLIC CYt. I NMICAl" REFLECTOR * ^AR* FIELD** V 

MOM/. MO.0/. .?{Si//°" " m "" '*““*”*•///• MOAO 

INPUT PARAMETERS- , 

WAVELENGTH OF ELECTRIC FIELD , „ t/ 

LOCATION OF COORDINATE ORIGIN WRT FEED (X,‘ Y,‘ 2) 3FB 3 

FEED ROTATION ANGLES (ALPHA , BETA, GAMMA ) .. ... ' • 3F8 3 

FOCAL LENGTH OF PARABOLIC CYLINDER.;;;:..: • F8 3/ 

ANGLE OF ROTATION ABOUT X-AXIS (PSI) .* IfB* 3) 


format ( 

I 

/• 

/• 

/• 

/• 

) 

FORMAT ( 


***** PROGRAM IN SINGLE PANEL MODE 
MINIMIM-Y POINT ON THE REFLECTOR 
MAXIM IM-Y POINT ON THE REFLECTOR 
MINIMtW-Z POINT ON THE REFLECTOR 
MAXIMUM-Z POINT ON THE REFLECTOR 


(XpY,Z) . \3P8.3 

(XpY,Z) • , 3P8. 3 

^fY,Z) 1 , 3F8. 3 

(XfY # Z) 9 , 3F8. 3 


/ 


APERTURE PLANE LOCATION (XC) • r7 

FEED SHADOW CENTER COORDINATES IN APERTURe'pL.* !2F7. 2 

HALF MAJOR AXIS OF FEED SHADOW • F7 j/ 

HALF MINOR AXIS OF FEED SHADOW *F7 3/ 

NUMBER OF PANELS IN REFLECTOR •xd/\ 

f^a?^;„^a t Io er pl;ne f - e ^ in one deg ***«“« 

FORMAT (' PHI-180 PLANE -•/) 

FORMAT (2X, SF14.fi) 

P I *DARCOS (-1 . ODO) 

PI 2»PI 4-PI 
PID2-0. 5*PI 
DTOR-PI/180. 

RTOD-180./PI 

RETURN 

END 


SUBROUTINE APRTUR (P, ICALL) 

IMPLICIT REAL*8 (A-H.O-Z) 

REAL*8 NHAT ( 3 ) ,NMAG 

INTEGER SURFCE 

COMMON /APR PRM/N PTPPL ,NPER I M 

COMMON /CONTRL/NOPT (3) ,NLIST,IOPT, ILIST(IOO) 

COMMON/FEED/EP (9 1 ) , ET (91 ) ,NP, NT 
COMMON /MATH/P I , PI 2, PID2, DTOR , RTOD 

COMMON/PARAMS/AORORF, BELLP, PSI , PLNPNT (3) , PLNORM ( 3 ) , FEED (3) . 

ALPHA, BETA, GAMMA, XCAM , XC , S URrCE, NPNL, NPOINT 


m 


DIMENSION A(3,3) ,AINV(3,3) ,B(3,2) ,BB(3,2) ,C(3) ,EI (3) ,ER(3) , 

P (5,2750) , SR (3) ,X (3) 

IF (ICALL.OT.l) GO TO SO 
ALPHAR-ALPHA*DTOR 
BETAR-BETA*DTOR 
GAMMAR-GAMMA*DTOR 

A (1 , 1 ) “DC OS (ALPHAR) *DCOS (GAMMAR) -DSIN (ALPHAR) *D8IN (BETAR) * 

. DSIN (GAMMAR) 

A ( 1 , 2) -DSIN (ALPHAR) *DCOS (GAMMAR) +DCOS (ALPHAR) *DSIN (BETAR) * 

. DSIN (GAMMAR) 

A ( 1 , 3) “-DCOS (BETAR) *DSIN (GAMMAR) 

A ( 2 , 1 ) — DS I N (ALPHAR ) *DCOS (BETAR ) 

A (2, 2)- DCOS (ALPHAR) *DCOS (BETAR) 

A (2, 3)- DSIN (BETAR) 

A ( 3 , 1 ) -DCOS (ALPHAR) *DSIN (GAMMAR) +DSIN (ALPHAR) *DSIN (BETAR) • 

. DCOS (GAMMAR) 

A ( 3 , 2) -DSIN (ALPHAR) *DSIN (GAMMAR) -DCOS (ALPHAR) *DSIN (BETAR) * 

. DCOS (GAMMAR) 

A (3, 3)- DCOS (BETAR) *DCOS (GAMMAR) 

DO 40 1-1,3 

DO 40 J-l, 3 
40 AINV (I , J) -A (J,I) 

NPERIM-4 

NPTPPL-2250 

50 IF (IOPT.EQ.O) CALL APRIN (P, ICALL) 

TMAX-0.0D0 

TMIN-PI 

PMIN-PI+PID2 

PMAX-PID2 

DO 65 I-1,NPERIM 

DO 60 J-l, 3 

60 X (J)-AINV(J, 1 ) *P (1, IJ+AINV (J , 2) *P (2, I) +AINV (J , 3) *P (3, 1) 

R-DSQRT ( (X (1 ) +FEED (1 ) ) **2+ (X (2) +FEED (2) ) **24- (X (3) 4-FEED (3) )**2) 
P (1 , I ) -DARCOS ( (X (3) +FEED (3) ) /R) 

SINTHT-DSIN (P(1,I)) 

IF (SINTHT.LT.l.D-10) SINTHT-1 . D-10 
P (2,1) -PI -DARSIN ( (X (2) 4-FEED (2 ) ) / (R*SINTHT) ) 

IF (P (1,1) .GT.TMAX) TMAX-P (1,1) 

IF (P (1,1) .LT.TMIN) TMIN-P (1,1) 

IF (P (2 , I ) . GT. PMAX) PMAX-P (2,1) 

IF (P (2, 1) .LT.PMIN) PMIN-P (2, 1) 

65 CONTINUE 

DELP-PMAX-PMIN 

DELT-TMAX-TMIN 

NP-DSQRT (DELP*DFLOAT (NPTPPL) /DELT) +1 . 0 
NP=( (NP-l)/2)*2+l 
ANGINC-DELP/ (DFLOAT (NP) -2 . 6) 

NTD2-DELT/ (2 . 0+ANGINC) +1.0 

NT=2*NTD2+l 

PMIN-PMIN-O. 8+ANGINC 

PMAX-PMAX+O . 8+ANGINC 

TCT= (TMAX+TMIN) /2. 0 

TM I N-TCT -DFLOAT (NTD2) +ANGINC 

TMAX-TCT+DFLOAT (NTD2) +ANGINC 

DO 95 J-l, NT 

DO 95 K-l.NP 

P (1 ,NPERIM+ (J-l) +NP+K)«TMIN+ (J -1 > +ANCINC 
95 P (2 , NPERIM+ (J -1 ) +NP+K) -PMIN+ (K-l ) +ANGINC 
NTNP-NT+NP 
NPOINT-NPER IM+NTNP 
TMIN-TMIN+RTOD 
TMAX-TMAX*RTOD 
PMIN-PMIN*RTOD 
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PMAX-PMAX*RTOD 

ANGINC-ANGINC*RTOD 

IP (10(1, ICALL) . BQ. 1 ) PRINT 107, TMIN,TMAX, PMIN, PMAX, 

. ANGINC,NTNP,NPOINT 

107 PORMAT(//' ILLUMINATION DATA-'// 

. ' THETA ILLUMINATION PROM '^9.3,' TO 

F9.3/ 

. ' PHI ILLIHINATION FROM ',F9.3,' TO 

P9. 3/ 

. ' INCREMENTAL ANGLE (DEG) ' ,F7. 4/ 

. 1 THEREFORE TOTAL NIMBER OF GENERATED RATS ',17 / 

. ' TOTAL NUMBER OF APERTURE PLANE POINTS ',17 ) 

IF (SURFCE.LT. 5) GO TO 114 
SNPSI-DSIN (PSI *DTOR) 

CSPSI-DCOS (PSI *DTOR) 

114 CALL FILL (P,NPOINT) 

DO 600 I«1 ,NPOINT 

SINP-DSIN (P (2, I) ) 

COSP-DCOS (P (2, 1) ) 

SINT“DSIN (P (1 , 1) ) 

COST"DCOS (P (1 , 1 ) ) 

BB(l,l)-SINT*COSP 
BB (2 , 1 ) »SINT*SINP 
BB (3, l)»COST 
BB ( 1 , 2) •+FEED (1 ) 

BB (2, 2) b+FEED (2) 

BB(3,2)-+FEED(3) 

CALL MULT32 (B,A,BB) 

GO TO (120,130, 140, 150, 160), SURFCE 
120 AR-0.0 

BRbB (1 , 1 ) *PLNORM (1 ) +3 (2 , 1 ) *PLNORM (2) +B (3, 1 ) *PLNORM (3) 

CR— (B ( 1 , 2) +PLNPNT (1 ) ) *PLNORM (1 ) 

. - (B (2, 2) +PLNPNT (2) ) *PLNORM (2) 

. - (B (3 , 2) +PLNPNT (3) ) *PLNORM (3) 

GO TO 180 

130 AR»B (1,1) **2/AORORF**2+ (B (2 , 1) **2*B (3 , 1) **2) /BELLP**2 

BRb- 2.0* (B (1, 1) *B (1 , 2) /AORORF**2+ (B (2, 1 ) *B (2, 2) +B (3, l) *B (3, 2) ) / 

. BELLP**2) 

CRbB(1,2)»*2/AORORF**2+(B(2,2)**2+B(3,2)**2)/BELLP**2-1.0 
GO TO 180 

140 AR-B(1,1)*B(1, 1)+B(2,1)*B(2,1)+B(3,1)*B(3,1) 

BRb- 2 . * (B ( 1 , 1 ) *B (1 , 2)+B (2, 1) *B (2 , 2) +B (3 , 1) *B (3 , 2) ) 
CR-B<l,2)*B(l,2)+B(2,2)*B(2,2)+B(3,2)*B(3,2)-AORORF*AORORF 
GO TO 180 

150 AR»S(2,1)*B(2,1)+B(3,1)*B(3,1) 

BR«-2. 0* (B (2, 1 ) *B (2, 2)+B (3, 1) *B(3, 2) +2. 0*AORORF*B (1, 1) ) 

CR*B (2, 2) *B (2 , 2)+B (3, 2) *B (3, 2)+4 . 0*AORORF*B (1 , 2) -4 .0*AORORF**2 
GO TO 180 

160 AR*B (3, 1) *B (3, 1) *CSPSI *CSPSI -2. 0*B (2 , 1) *B (3, 1) *CSPSI*SNPSI 
. +B (2 , 1 ) *B (2 , 1) *SNPSI *SNPSI 

BR— 2.0* (B (3, 1 ) *B(3, 2) *CSPSI*CSPRI 
. - (B(3, 1) *B(2, 2) +8 (2, 1) *B(3,2) ) *CRPSI *3NPSI 

. +B(2,1)*B(2,2) *SNPSI ‘SNPSI+2. 0*AORORF*B (1,1) ) 

CR-B(3,2)*B(3,2)*CSPSI*CSPSI-2.0*B(2,2)*B(3,2)*CSPSI*SNPSI 
. +B (2 , 2) *R (2 , 2) *SNPS I *SNPSI+4 . 0* AORORF* (B ( 1 , 2 ) -AORORF) 

180 IF (AR.LT.1.0D-10) R— CR/BR 
IF (AR.LT.1.0D-10) GO TO 189 
R« (-BR+DSQRT (BR*BR-4 . *AR*CR) ) / (AR+AR) 

189 CONTINUE 

X0«B (1 , 1) *R-B (1 , 2) 

Y0«B (2 , 1 ) *R-B (2 , 2) 

Z0“B (3, 1) *R-B (3, 2) 

219 GO TO (220, 230, 240, 250, 260), SURFCE 
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220 NHAT (1 ) -PLNORM ( 1 ) 

NHAT (2) -PLNORM (2) 

NHAT ( 3 ) -PLNORM ( 3 ) 

GO TO 288 

230 NHAT(l)— XO*BBLLP**2/DSQRT(XO**2*BELLP**4+<YO**2+ZO**2)*AORORP**4) 

NHAT (2 ) * -Y 0 *AORORF* * 2 /DSQRT (X0**2*BELLP**4+ (Y0**2+Z0**2) * 

. A0R0RF**4) 

NHAT (3) *-ZO*AORORF**2/DSQRT (X0**2*BELLP**4+ (Y0**2+Z0**2) * 

. A0R0RF**4 ) 

GO TO 288 

240 NHAT (1 ) --XO/AORORF 
NHAT (2 ) — YO/AORORP 
NHAT (3 ) --ZO/AORORF 
GO TO 288 

250 NHAT (1 ) -2 . 0*AORORF/DSQRT (4 . 0*AORORF**2+YO**2+ZO**2) 

NHAT (2 ) - -YO/DSQRT (4 . 0*AORORF**2+YO**2+ZO**2) 

NHAT (3 ) • -ZO/DSQRT (4. 0*AORORF**2+YO**2+ZO**2) 

GO TO 288 

260 NMAG-DSQRT(4.0»AORORF*AORDRFMZO*CSPSI*.,NPSI-YO*SNPSI*SNPSI)**2 
. + (Y 0 *SNPS I *CSPS I -Z 0 *CSPS I *CSPS I ) * *2 ) 

NHAT (1 ) - 2 . 0*AORORF/NMAG 
NHAT(2)-SNPSI* (ZO*CSPSI-YO*SNPSI)/NMAG 
NHAT (3 ) »CSPSI * (Y 0 *SNPS I -Z 0 *CSPS I ) /NMAG 

288 SCALAR-2. 0* (B (1 , 1) *NHAT (1 )+B (2, 1) *NHAT (2) +B (3, 1) *NHAT (3) ) 

DO 295 L-1,3 

295 SR (L)*B (L, 1) -SCALAR*NHAT (L) 

ETI -P (3 , 1 ) /R 
EPI -P (4 , 1 ) /R 

C(1)-C0ST*C0SP*ETI-SINP*EPI 

C(2)»C0ST*SINP*ETI-K:0SP*EPI 

C(3)— SINT*ETI 

DO 400 N-1,3 

El (N) -0.0 

DO 400 M-1,3 

400 El (N)-EI (N)+A(N,M) *C(M) 

SCALAR-2.0* (El (1 ) *NHAT (1 ) +EI (2) *NHAT (2)+EI (3)*NHAT(3)) 

DO 500 K-1,3 

500 ER (K ) -SCALAR*NHAT (K ) -El (K) 

Y-Y0+ (XC-XO) *SR (2) /SR (1 ) 

7-20+ (XC-XO ) *SR (3 ) /SR (1 ) 

D-DSQRT ( (XC-XO) * (XC-XO) + (Y-YO) * (Y-YO) + (Z-ZO) * (Z-ZO) ) 

PHASE-PI 2* (R+D) /XLAM+P (5,1) 

P (1 , I ) -Y 
P (2 , 1 ) -Z 
P(3, 1 ) — ER ( 2 ) 

P (4 , 1 ) -ER ( 3 ) 

P (5, D-PHASE 

600 CONTINUE 
RETURN 
END 


SUBROUTINE APRIN (P, ICALL) 

IMPLICIT REAL*R (A-H.O-Z) 

INTEGER SURPCE 
COMMON/APRPRM/NPTPPL, NPERIM 
COMMON/CONTRL/NOPT (3) ,NLIST, IOPT, ILIST (100) 

COMMON/PARAMS/AORORP, BE LLP, PS I , PLNPNT (3 ) , PLNORN ( 3 ) ,FEED(3) , 
ALPHA, BETA, GAMMA, XLAM,XC, SURPCE, NPNL.NPOINT 
DIMENSION P (5 , 2750 ) 

READ (1 , 10) NPERIM, SURFCE,NPTPPL 
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10 FORMAT (31 5) 

IF (NPERIM. LE. 2) CO TO 2S0 

IF (NPERIM. OT. 40) GO TO 260 

IF (SURFCE.Cr.5) GO TO 250 

IF (NPTPPL.GT.2S00) GO TO 270 

IF ( (NPERIM*S URFCE) .LB.O) GO TO 250 

READ(1,20) ((P(I,J> ,1*1 >3) »J*1 , NPERIM) 

20 FORMAT (3F10. 4) 

28 GOTO (30,40,50,50,60) ,8 URFCB 

30 READ(1 , 35) PLNPNT,PLNORM 

35 FORMAT (6F10.4) 

GO TO 100 

40 READ (1,4 5) A0R0RP,8BLLP 

45 FORMAT (2F10.4) 

GO TO 100 

50 READ (1 , 55) AORORF 

55 FORMAT (F10. 4) 

GO TO 100 

60 READ(1,65) AORORF, PSI 

65 FORMAT (2F10. 4) 

100 CONTINUE 

199 IF (10(1, ICALL) . EQ. 0) RETURN 
PRINT 200, ICALL 

200 FORMAT ( ' 1 ' , 35X, 'REFLECTOR PANEL NUMBER', 14) 

GO TO (320,330, 340, 350, 360) .SURFCE 

250 PRINT 252, ICALL 

252 FORMAT (///' ********** INPUT ERROR ON CARD ONE FOR PANEL NUMBER', 

. 14,' EXECUTION TERMINATING ***•**•***•) 

STOP 

260 PRINT 262, ICALL 

262 FORMAT (///• ********** STORAGE DOES NOT EXIST FOR NUMBER OP', 

. * PERIMETER POINTS SPECIFIED - PANEL', 14,' *•****•*•*•) 

STOP 

270 PRINT 272 ICALL 

272 FORMAT (///' ********** MAXIMUM ILLUMINATION REQUEST IS 2500', 

. ' RAYS - PANEL', 14,' **********•////) 

NPTPPL-2500 
GO TO 28 

320 PRINT 401, PLNPNT, PLNORM,NPERIM 
RET URN 

330 PRINT 402, AORORF, BELLP,NPERIM 

RETURN 

340 PRINT 403, AORORF, NPERIM 

RETURN 

350 PRINT 404, AORORF, NPERIM 

RETURN 

360 PRINT 405, AORORF, PSI, NPERIM 

RETURN 

401 FORMAT (///10X, 'PANEL IS A PLANAR SURFACE',/// 

' A POINT ON THE REFLECTOR SURFACE (X,Y,E) *,3F7.2 

. /' COMPONENTS OF UNIT NORMAL TO SURFACE (X,Y,Z) *,3F7.2 

. /' NUMBER OF USER-SUPPLIED EDGE POINTS ',17) 

402 PORMAT(///l OX, 'PANEL IS AN ELLIPTICAL SECTION',/// 

. ' MAJOR AXIS OF ELLIPTICAL REFLECTOR *,F7.2/ 

. * MINOR AXIS OF ELLIPTICAL REFLECTOR ',F7.2/ 

. ' NUMBER OF USER-SUPPLIED EDGE POINTS ',17) 

403 FORMAT(///l OX, 'PANEL IS A SPHERICAL SECTION'/// 

' RADIUS OF REFLECTOR SPHERE ',F7.2/ 

. ' NUMBER OF USER-SUPPLIED EDGE POINTS ',17) 

404 FORMAT (///l OX, 'PANEL IS A PARABOLIC SECTION',/// 

• FOCAL LENGTH OF THE PARABOLA .....*,F7.2/ 

. ' NUMBER OF USER-SUPPLIED EDGE POINTS ',17) 

405 FORMAT (///l OX, 'PANEL IS SECTION OF A PARABOLIC CYLINDER'/// 
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• FOCAL LENGTH OF THE PARABOLA ' ,F8. 3/ 

• FOCAL LINK ROTATION FROM Y-AXIS (FS1) • ,F8. V 

• NUMBER OF USER-SUFFLIEO EDGE POINTS ',I7> 


SUBROUTINE FILL(P,NPT) 

IMPLICIT R£AL*8 (A-H,0-Z) 

COMMON/FEED/EP (91) ,ET(91) ,NP,NT 
COMMON/MATH/P I , PI 2 , PID2 , DTOR , RTOD 
DIMENSION P(S,NPT) 

DO 100 1*1, NPT 

PROJX*D8IN (P (1 , I)) *DCOS (P <2, 1 ) ) 

PROJY-O.ODO 

IP (DABS (P (2, 1) -PI). OT. 1.00-5) 

.PROJY-DCOS (DATAN (DCOS (P (1 , I) )/DSlN (P (1 , I) )/DSIN (P (2, 1) ) ) ) 
AN6LX-DARC0S (DABS (PROJX) ) *RTOD 
LO-ANGLX+l.ODO 
IHI»L04l 

PPFLD* (ANCLX-DFLOAT (LO-1) ) * (EP (IHI ) -EP (LO) ) 4EP (LO) 

TPFLD* (ANGLX-DFLOAT (LO-1 ) ) * (ET (IHI)-ET(LO) )4ET (LO) 

COSE 2»PR0JY *PRO JY 
S INE2»1 , 0D0-C0SE2 
P (3, I)»O.ODO 

P (4 , I) *PPFLD*TPFLD/ (DSQRT (TPFLD*TPFLD*C0SE2+ 

. PPF LD*PPF LD*S INE 2 ) ) 

P (5, 1)*O.ODO 
100 CONTINUE 
RETURN 
END 
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SUBROUTINE QUANTZ (P,NPERIN, ICALL) 
IMPLICIT REAL‘8 (A-H.O-Z) 

INTEGER SURFCE 


COMMON/BLOCKC/YCBL, ZCBL, HPMABL, HFHIBL 
COMMON/CONTRL/NOPT (1) .NLIST, I OPT, ILIST (100) 

COMMON /DIMENS/YDIM,ZDIM,YCT, ZCT 
C0MM0N/EXTENT/YMIN,YMAX,ZM1N,ZMAX 
COMMON/FEED/EP (91) ,ET(91) ,NP,NT 

COMMON/PARAMS/AORGRF, BELLP, PSI , PLNPNT (3) .PLNORM ( 3) ,FEED(3) , 

. ALPHA .BETA, GAMMA , X L AM , XC , S UR FC E , NPN L , NPO INT 

DIMENSION P (5, 2750) .PINT (5) , POLD (5) ,P8LK (5) ,PR5 (5,41) ,Z(2,101) 
NBARS-NP-2 
YMIN-1.0D410 
YMAX— 1.0D+10 
ZMIN-1.0D+10 
Z MAX— 1.00*10 


NOS*2*NBARS 

CALL SETM ( l . 0D+20,Z,NOSt 
DO 20 1*1 ,NPERIM 

IP (P (1,1) .CT. YMAX) YMAX*p (1,1) 
IP (P (1 , I ) . LT. YMIN) YMIN*P (1,1) 
IP (P(2,I) .GT.ZMAX) ZMAZ*P(2,t) 
IP (P(2,I) .LT.ZMIN) ZMIVap (2, I ) 
CALL MOVEM (P (1 , 1 ) , PR5 (1, 1 ) , 5) 
YDIM “YMAX -YMIN 
YCT- (YMAX+YMIN) /2. 
ZDIM-ZMAX-ZMIN 




fJ 0, 


0 ]{ 




ZCT* (ZMAX+ZMXN)/2. 

GRID* (YNAX-YRIN) / (DPLOAT (MBARB) -0.4) 
GRXDLO-YMIN+ORID/S.ODO 
GRIDHIaYRAX-CRID/S.ODO 
31 tB»*NPERIM+l 
NDEX-NPBRIN 
DO 100 I-IBON,NPOINT 

ip (Pd.n.or.YNAX) go to m 

IP (P(1,I) .LT.YHIN) GO TO 01 
NCR ID* (P (1,I)-GRIDLO)/GRID+O.S 
P (1 , 1) -GRIDLO+DPLOAT (MAID) *CRID 
CALL WOVEN (P (1, t) , P (1, IHIDEX) ,5) 

GO TO 100 
9S NDEX-NDSX+1 
100 CONTINUR 

NPOINT-NPOINT-NDEX 
CALL PT80RT (P, S,NPOIHT) 

IP (IOPT.BO.1) GO Tlf 422 

CALL WOVEN (PRS (1 , 1 ) , ; R5 (1 , NPERIN41 ) f 5) 

KDBX-2 

Y2*PRS (1*1) 

Z2*PR5 (2*13 
200 Yl-PR5(l,KOEX) 

Z1*PR5 (2,KDEX) 

IP (DABS (Y1-Y2) . LT, 1. OD-S) GO TO 400 
SL0PB*<Z1-Z2)/<Y1-Y2) 

R*Z2-SLOPB # Y2 
IP (Y1-Y2) 220,230,230 
220 /HI-Y2 
YLO-Yl 
GO TO 240 
230 YHI*Y1 
YLO-Y2 

240 INDEX* <YLO-CRIDLO)/GRID4l.O 
250 INDEX-INDBX41 

YQ*GRIDLO+DPLOAT (INDEX-1 ) ‘GRID 
IP (YO.OT.YHI) GO TO 400 

1 LOAD* 1 

ZEE*SLOPE*YQ+B 

IP (Z(l, INDEX) .LT.1.0D+10) I LOAD- 2 
7. ( I LOAD, INDEX )*Z EE 
GO TO 250 
400 Y2*Y1 

2 2*21 

KDEX-KDEX+1 

IF (KDEX. LE.NPERIN+1 ; GO TO 200 
DO 420 I-l.NBARS 

l' ( (2(1, 1)42 (2,1) ! CT.l. 00*10) GOTO 1005 
.2(2, I) -Z (1,1)1 410,420,420 
'•’0 >.2*2 (2,1) 

2(2,1) 'Z(l.I) 

Z ( 1 , 1 ! *Z 2 
420 CONTINUE 
n T’> 444 

422 HPMAEX*YDIN/2, 0D0 
H PMIEX-ZDIN/2.0D0 
DO 430 I-l , NBARS 
Y-GRIDL04DPL0AT(I-1)*GRID 
ZZ-HPNI EX*DSQRT ( 1 . ODD- ( (Y -YCT) /HPNAEX ) **2 ) 
2(1,1)*- ZZ+ZCT 
Z (2, 1)* ZZ+ZCT 
410 CONTINUE 
444 L-0 
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470 

480 


SOI 

S10 

SIS 

701 

704 


711 


720 

72S 

800 

900 


9 f >0 


N-i 

CALL SETM (0. 0, PBLK, 5) 

YQ-P (1*1) 

I DCX-DINT ( (YQ-GR IDLO) /GAID>1 .001) 

DO 900 I -1, NPOINT 
IP <P(1,I) .EQ.YQ) CO TO 481) 

IP (L.OT.2) GO TO 470 

N-N-L 

L-0 

YQ-P (1,1) 

tDEX-DINT ( (YO-CRIOLO) /CHID* 1.001) 

PBLK (1 )-P (1 , 1 ) 

P8LK (2) -P (2, 1 ) 

TEST— 1.0 

IF (P(2,I).E0.Z(1« IDEX) .OR.P (2, I ) . CO. 2(2. IDEX) ) TEST- 0.0 
IP (P (2, 1 ) .OT. 2(1, IDEX) .AND. P (2. I ) . LT.Z (2, IDEX) ) TEST-1. 0 
TESTS L-HFMABL*HFMABL*NPMIBL*HPMIBL 

-HFMABL*HFMABL* (P (2 , I ) -2CBL) • (P (2 , I) -2CBL) 
-HFHIBL*HPMIBL- (P (1 , I ) -YCBL) * (P (1 , I ) -YCBL) 

IP (TEST) 701,501,501 
IP (TESTBL.LE.O.O) CO TO 510 
CALL NOVEM (PBLK, P (1 , N) , S) 

CO TO 515 

CALL NOVEM (P(1,I),P(1,N),5) 

N-N + l 
L«L*1 

TP (TEST. EQ. 0.0) CO TO BOO 
IP (L.EQ.O) CO TO 800 

IP (TEST*TESTO) 704,800,800 
CALL INTPL (POLO, P (1,1), PINT, 2(1, IDEX)) 

NCHG-0 

IP (TEST. LT. 0.0) CO TO 711 
CALL NOVEM (P ( 1 , N — 1 ) ,P(1,N) ,5) 

NCHC-1 

CALL NOVEM (PINT, PBLK, 2) 

TESTBL-HPMABL*HFMABL*HPNIBL*HFNIBL 

. -HPMABL*HFMABL* (PINT (2 ) -2CBL) * (PINT (2) -2CBL) 

-HFMIBL'NFMIBL* PINT (1 ) -YCBL ) * (PINT (1 ) -YCBL) 
IP (TESTBL.LE.O.O) CO TO 720 
CALL NOVEM (PBLK, P (1 , N-NCHG) , 5) 

CO TO 725 

CALL NOVEM (PINT , P (1 ,N NCHG) , 5) 

N«N+1 
L-L + l 

CALL MOVEM (P (1 , I) ,POLD, 5) 

TESTG-TEST 
CONTINUE 
NPOINT-N -1 
NLOC-S*NPERIM 

CALL MOVEM(PR5,P(l.N),NLOC; 

IF (IO(4,l).EQ. 1) WRITE (",<»•> ) NPOINT,ICALL 

IF (10(4, 1) .EO. 1) WRITE (7 9-2) ( (P (l , J ) ,1-1,5) , J-l , NPOINT) 

IP (10(1 , ICAt.L) .E0.0) RETURN 

PR I NT 950, YMIN, YMAX, 2NIN, ZMAX, GRIDLO, GRIDNI , CR l D.NBARS, NPOINT 


Format (// ■ quantizing data-*// 

' POINT PATTERN EXTENTS ON APERTURE PLANE YNIN-*,P7.2/ 

YNAX-* ,P7. 2/ 

' ZNIN-'.PT.I/ 

’ ZNAX-* ,P7. 2/ 

' GRID RANGES PROM *,P8.3.' TO 

ra.v 

• SPACING BETWEEN GRID BARS IS ,.*,P*.4/ 

' THEREPORE NUMBER OP GRID BARS ',!*/ 



. * NUMB KM OF FOINTB SUPPLIED TO RADFAT * # t? ) 

951 FORMAT (311) 

952 FORMAT (5016. 10) 

RBTURN 

1005 •'HINT 1006 

1006 FORMAT (/// 

. «••••••••* : ARRAY FOR FOINT RKBIDKNCC NOT PILLED CORRECTLY 

. •••••••••»', //10X,'- STOF EXECUTION -•) 

STOP 

END 


SUBROUTINE I NT CRT (P, MAJOR, AMAJOR , MINOR , FI CLOY , FIELD! ) 
IMPLICIT REAL*! (A-M,0-Z) 

REAL*! MAJOR, MINOR 

COMPLEX*16 CTEMP.CZ 1,CZ2,CY1,CY2,TSZ ,TSY, DZI , DYI.ZIOLD, YIOLD, 
. Z I, YI.FLDZ, PLOY, FIELDS (200) ,FIELDY(200) 

INTEGER SURFCE 

COMMON /MATN/P I , PI 2, PID2,DTOR, RTOD 

COMMON/PARAMS/AORORF, BELLP, PR I , PLNPNT ( J ) , PLNORM (1 ) , FEED ( 3 ) , 
ALPHA , BETA , GAMMA , XLAM , XC , 8 URFCE, NPN L , NPOINT 
DIMENSION AMINOR (3) ,P (5, 2750) 

DATA HPHI , HTHTA/5HPHI , 5HTNETA/ 

SEN-999.0 
NPARTS-7 
R PART- 1 ./N PARTS 
ZLAM-PI 2/XLAM 

CALL SETH (SEN, P (1 , NPOINT+1 ) , 5 ) 

DEC- AMAJOR 

DEGR-DEC*DTOR 

DLOR -AMINOR ( 1 ) *DTOR 

DICR-AMINOR (I) *DTOR 

DSTOPR -AMINOR (2) *DTOR*DICR«0 . 5 

NTH-0 

D-DLOR 

IF (MAJOR. NE. HPHI) GO TO 1400 
400 COSP-DCOS (DEGR) 

SINP-DSIN (DEGR) 

COST-DCOS (D) 

SINT-DSIN (D) 

GO TO 3425 
>400 COSP-DCOS (D) 

SINP-DSIN (D) 

COST-DCOS (DEGR) 

JINT-DSIN (DEGR) 

1425 NTH-NTHM 

CTSP-COST*SlNP 
ZK-Z LAM*COST 
YK-Z LAM*SINP*5 INT 
(OLD- l 
INEW-2 

PLDY-(0. 0,0.0) 

PLDZ-(0. 0,0.0) 

Y OLD-SEN 
YI- (0.0,0. 0) 

Zl-(0. 0,0.0) 

1450 CONTINUE 

IP (P (1 , IOLO) .NE. P(l.tNEW)) CO TO 4000 
Z-P (2, IOLD) 

ERY-P (3, IOLD) 




PH-P(5,I0LD) 

D2« (P (2, INEW) ~Z) *RPART 
DEI\/~ (P (3, INEW) -ERY) *RPART 
DER2* (P (4 , TNEW) -ERZ) *RPART 
DPH* (P (5, INEW) -PH) *RPART 
CTEMP«CDEXP (DCMPLX (0 . 0D0 , ZK*2-PH ) ) 
CZ l"ERZ*COSP # CTBHP 
CY1» (ERY*SINT+ERZ*CTSP) *CTBMP 
TSY* (0 • 0, 0. 0) 

TSZ*{0. 0,0.0) 

DO 3700 N*1 , NPARTS 
Z*Z+DZ 

ERY*RRY+DERY 

ERZ-ERZ+DERZ 

PH*PH+DPH 

CTEMP-CDEXP (DCMPLX (0 . 000, ZK*Z-PH) ) 
CZ2»ERZ*COSP*CTEMP 
CY2* (ERY*SINT+ERZ*CTSP) *CTEMP 
TSZ*TS*+CZ1+CZ2 
TSY*TSY+CY1+CY2 
CZ 1*CZ 2 
CY1*CY2 
3700 CONTINUE 

Z I *Z I +TSZ * ( . 5*DZ ) 

YI-YI+TSY* ( . 5*DZ ) 

3900 IOLD=IOLD+l 
INEW=INEW+1 
GO TO 3450 
4000 CONTINUE 

YNEW=P (1 , IOLD) 

IF (YOLD. EQ. SEN) GO TO 4400 
4200 DZI* (Z I -Z IOLD) *RPART 
DYI* (YI-YIOLD) *RPART 
DY= (YNEW-YOLD) *RPART 
CTEMPCDEXP (DCMPLX (0. 0D0, YK*YOLD) ) 
CZ1=ZI0LD*CTEMP 
CY 1 *Y I0LD*CTEMP 
TSY* (0.0, 0.0) 

TSZ*(0. 0,0.0) 

DO 4300 N=1 , NPARTS 
YOLD*YOLD+DY 
Z IOLD*Z IOLD+DZ I 
YI0LD*YI0LD+DYI 

CTEMP*CDEXP (DCMPLX (0 . ODO, YK*Y0LD) ) 
CZ2=ZIOLD*CTEMP 
CY 2=Y IOLD*CTEMP 
TSZ=TSZ+CZ 1+CZ 2 
TSY*TSY+CY1+CY2 
CZ1-CZ2 
CY 1*CY 2 
4300 CONTINUE 

F LDZ =FLDZ+TSZ * (. 5*DY) 
FLDY*FLDY+TSY* (. 5*DY) 

4400 CONTINUE 
YOLD“YNEW 
ZIOLD-ZI 
Y I QLD*Y I 
YI*(0. 0,0.0) 

ZI=(0. 0,0.0) 

IF ( P ( l , I N EW ) . NE • SEN) GO TO 3900 

FIELDY (NTH)-FLDY 
FIELDZ (NTH)-FLDZ 
D=D+DICR 
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IF (D.GT.DSTOPR) GO TO 5000 
IF (MAJOR. EQ.HPHI) GO TO 400 
GO TO 3400 
5000 CONTINUE 
RETURN 
END 
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• OUTPUT FOR THE EXAMPLE 
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